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Preface 

Although  this  study  was  begun  with  the  thought  in  mind  of 
providing  an  unbiased  evaluation  of  various  describing-function 
generating  schemes  (unbiased,  since  I  had  no  hand  in  any  of  them) , 
it  has  ended  up  as  an  introduction,  to  a  describing-function  genera¬ 
ting  scheme  of  my  own,  as  well.  I  should  not  really  call  the  correc¬ 
ted-  conventional  describing  function  my  own,  since  it  is  really  only 
an  unconventional,  "corrected"  form  of  the  conventional  describing 
function  developed  by  Kochenburger  and  others,  and  its  only  theore¬ 
tical  justification  is  that  it  is  based  on  the  conventional  tech¬ 
nique.  Let  me  assure  you,  however,  that  the  basic  linear  functions, 
used  as  being  representative  for  the  purpose  of  the  study,  were 
chosen  before  the  idea  of  correcting  the  conventional  describing 
function  came  to  me,  and  therefore  they  in  no  way  (intentionally, 
at  least)  were  chosen  because  they  favor  the  corrected-conventional 
technique.  As  a  matter  of  fact,  I  first  realized  the  possible  ben¬ 
efit  of  such  a  correction  while  trying  to  improve  on  another  uncon¬ 
ventional  technique  presented  by  Gibson  called  the  new  rms  describ¬ 
ing  function.  And  it  was  on  his  examples  that  I  first  had  success. 

I  wish  to  express  my  appreciation  for  the  encouragement,  guid¬ 
ance,  and  assistance  that  my  Faculty  Thesis  Advisor,  Professor  J.  J. 
D'Azzo  provided  me,  as  well  as  for  the  original  topic.  I,  further¬ 
more,  would  like  to  thank  Hr.  R.  0.  Anderson  of  the  Flight  Control 
Division  (RAID)  for  sponsoring  my  work  and  also  providing  guidance 
and  advice.  Some  of  the  greatest  tangible  assistance,  though,  came 
from  Mr.  H.  E.  Petersen  of  the  Analysis  Branch  of  the  Digital  Com- 
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puter  Division  (R&TD)  vho  programmed  and  ran  the  describing  function 
equations  on  the  digital  computer  so  that  they  could  be  plotted,  and 
Mr.  David  K.  Bowser  of  the  Control  Criteria  Branoh  of  the  Flight 
Dynamics  Laboratory,  Flight  Control  Division  (B&TD),  vho  envisioned, 
programmed,  and  ran  the  digital  computer  verification  of  the  analog 
computer  results.  Last,  but  not  least,  I  thank  ray  wife  for  typing 
this  thesis. 


Robert  R.  Rankine,  Jr, 
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Abstract 

A  describing  function  is  an  amplitude-dependant  generalized 
transfer  function  of  a  nonlinearity  which  can  be  used  to  repre¬ 
sent  the  nonlinearity  when  its  input  is  approximately  sinusoidal. 

The  simpler  and  more  prominent  or  accurate  describing-function 
generating  schemes  include  (l)  the  conventional,  (2)  the  minimum 
average  error,  (3)  the  equivalent  gain,  (4)  the  new  rms,  and  (5) 
the  corrected-conventional.  An  evaluation  of  these  various  schemes, 
based  on  the  accuracy  with  which  the  describing  functions  they  pro¬ 
duce  can  predict  the  amplitude  and  frequency  of  self-sustained 
oscillations  in  a  nonlinear  system,  reveals  that  the  corrected- 
conventional  describing  functions  are  more  accurate  for  single¬ 
valued  nonlinearities. 
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AN  EVALUATION  OF  SELECTED  DESCRIBING  JUNCTIONS 
OF  CONTROL  SYSTIM  NONLINEARITIES 


I.  Introduction 


A  "generalized  transfer  function"  is  an  approximately  equiv¬ 
alent  linear  transfer  function  of  a  nonlinearity.  It  is  used  to 
represent  the  nonlinearity  in  control  Bystem  analysis,  but  is  appli¬ 
cable  only  when  the  input  to  the  nonlinearity  is  sinusoidal.  If  the 
nonlinearity  is  independent  of  the  inpit  frequency, <*J  ,  and  is  solely 
dependent  on  the  input  amplitude,  X,  the  generalized  transfer  func¬ 
tion  is  known  as  a  "describing  function"  or  moire  specifically,  a 
"sinusoidal  describing  function"  (Ref  4:405)*  One  derives  the  de¬ 
scribing  function  of  a  particular  nonlinearity  by  using  some  scheme 
to  choose  a  sinusoid  that  will  closely  represent  the  output  of  the 
nonlinearity  in  some  sense  (Ref  7:153).  The  describing  function  then 
becomes  the  ratio  of  the  amplitude  and  phase  of  the  chosen  output 
sinusoid  to  the  amplitude  and  phase  of  the  input  sinusoid;  that  is, 


Describing  Function 


N 


Amplitude  and  Phase  of  the  Equivalent 
Sinusoidal  Output  of  the  Nonlinearity 
Amplitude  and  Phase  of  the  Sinusoidal 
Input  to  the  Nonlinearity. 


The  foregoing  definition  of  "describing  function"  is  far  more 


general  than  that  usually  encountered  in  literature  on  the  subject. 
In  the  past,  definitions  of  the  function  have  usually  encompassed 
only  the  describing  function  scheme  developed  by  Kochenburger 
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(Ref  8:270ff)  and  others^  based  on  the  Fourier  Series  representation 
of  the  output  of  the  nonlinearity.  Many  other  schemes  for  repre¬ 
senting  nonlinearities  have  been  developed  since  Kochenburger' s  work, 
however,  and  limiting  the  definition  to  this  "conventional”  describ¬ 
ing  function  is  apt  to  cause  confusion  when  one  tries  to  define  the 
others.  Thus,  the  term  "describing  function"  will  be  used  in  its 
most  general  sense  in  this  paper,  and  the  describing  function  orig¬ 
inated  by  Kochenburger,  which  represents  the  output  of  the  nonlin¬ 
earity  by  the  first  harmonic  of  the  Fourier  Series  expansion  of  the 
output  waveform,  will  be  referred  to  as  the  "conventional  describ¬ 
ing  function." 

Since  there  are  many  schemes  for  producing  describing  functions, 
the  question  arises:  "Which,  if  any,  is  best?"  The  purpose  of  this 
paper  is  to  answer  that  question. 

The  Problem 

The  problem  is  to  evaluate  comprehensively  the  accuracy  with 
which  the  various  describing  functions  can  predict  the  amplitude  and 
frequency  of  self-sustained  oscillations  in  a  nonlinear  feedback  con¬ 
trol  system.  When  this  is  done,  it  may  be  possible  to  choose  a 
describing-function  generating  scheme  which  produces  describing  func¬ 
tions  superior  to  all  of  the  others,  or  it  may  be  found  that  each  type 
of  nonlinearity  favors  a  different  scheme.  In  either  case,  the  con- 


Kochenburger  in  the  United  States,  Dutilh  in  France,  and  Goldfarb  in 
Russia  are  generally  given  equal  credit  for  having  discovered  the 
method  since  their  works  were  made  known  about  the  same  tine. 
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trol  system  engineer  shall  be  able  to  obtain  greater  accuracy  in 
predicting  the  efiecxs  of  a  nonlinearity  in  a  feedback  control  sys¬ 
tem  as  a  result  oi  this  study.  This  paper  is  not  intended  to  be  a 
primer  on  the  describing  function  method  of  analysis,  but,  rather,  a 
comprehensive  evaluation  of  a  fev  of  the  describing  functions  which 
may  be  used  in  that  analysis.  Therefore,  if  the  reader  is  unfamil¬ 
iar  with  the  general  describing  function  technique  for  stability  and 
limit  cycle  determination  by  direct  polar  (i.e.,  Hyquist)  plots,  it 
is  recommended  that  he  first  consult  the  references  (e.g. ,  Ref  1:442- 
444). 

Describing  Function  Generating  Schemes  to  be  Studied.  Since  it 
would  be  difficult,  if  not  impossible,  to  research  and  include  every 
scheme  for  generating  describing  functions  ever  invented,  only  the 
more  simple  and  prominent  or  accurate  methods  will  be  included.  These 
are  (l)  the  conventional  describing  function  originated  by  Kochen- 
burger  which  expands  the  output  of  the  nonlinearity  in  a  Fourier 
series  and  uses  the  first  Fourier  harmonic  as  the  equivalent  sinusoi¬ 
dal  output  (Rei  6:270ff)>  (2)  the  minimum  average  error  describing 
function  originated  by  Gibson  which,  as  its  name  indicates,  mini¬ 
mises  the  average  error  between  the  equivalent  sinusoid  and  the  out¬ 
put  ox  the  nonlinearity  (Ref  2:3&l);  (3)  the  equivalent  gain  origi¬ 
nated  by  r'rince  which  uses  the  maximum  amplitude  of  the  output  of  the 
nonlinearity  as  the  amplitude  of  the  equivalent  sinusoid  (Ref  10:217); 
(4)  she  new  ms  describing  function  originated  by  Gibson  which  is 
base!  on  equating  the  rms  value  of  the  actual  nonsinusoidal  output  of 
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the  nonlinearity  to  the  rms  value  of  the  equivalent  sinusoidal  out¬ 
put  (Ref  3:1321) j  and  (5)  the  correoted-conventional  describing  func¬ 
tion  developed  "by  the  author  which  includes  the  effects  of  Fourier 
harmonics  higher  than  the  first  by  taking  as  the  amplitude  of  the 
equivalent  output  sine  wave,  the  square  root  of  the  sul  of  the  squares 
of  the  amplitudes  of  the  first  and  third  Fourier  harmonics. 

One  of  the  more  prominent  methods  of  improving  the  accuracy  of 
the  conventional  describing  function,  the  use  of  Johnson's  correction 
terms  (Ref  6:l69ff)>  will  not  be  included.  The  calculation  of  the 
amplitude  correction  term  is  too  complex  to  be  very  practical  to  the 
engineer.  One  important  result  obtained  by  Johnson,  however,  will  be 
understood  throughout  this  study:  that  the  first  correction  term  for 
the  fundamental  frequency  is  zero,  and  that  therefore  the  describing 
function  frequency  prediction  is  generally  quite  accurate.  Because 
of  this  fact,  this  study  will  be  primarily  aimed  at  improving  the 
accuracy  of  the  amplitude  prediction;  nevertheless,  variations  in 
the  predicted  frequencies  will  result  when  hysteresis  is  studied  by 
the  different  schemes. 

Ronlinearities  to  be  Studied.  In  order  to  study  how  well  a 
particular  describing-function  generating  scheme  can  produce  a 
describing  function  which  accurately  predicts  a  system  oscillation, 
it  is  necessary  to  choose  a  number  of  nonlinearities  to  which  to 
apply  the  scheme.  The  most  often  encountered  nonlinearities  are  con¬ 
sidered  the  most  important,  and  the  following  nonlinearities  are,  in 
the  opinion  of  the  author,  the  most  often  encountered:  (  3  >  dead 
zone,  (2)  saturation,  (3)  dead  zone  combined  with  saturation,  (4)  the 
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ideal  relay,  (5)  the  relay  combined  with  hysteresis,  (6)  the  relay 

combined  with  dead  zone,  and  (7)  the  relay  combined  with  both  dead 

zone  and  hysteresis.  Each  of  these  nonlinearities  has  an  output 

which  is  both  odd  periodic  and  odd  harmonic  for  sinusoidal  inputs, 

o 

or  at  least  can  be  considered  so  by  phase-shifting  the  axes.  Each 
of  the  5  describing-function  generating  schemes  chosen  will  be 
applied  to  each  of  the  above  7  nonlinearities  producing  35  different 
describing  functions  to  be  studied. 

Linear  Systems  to  be  Studied.  It  is  also  necessary  to  choose 
some  arbitrary,  representative  linear  feedback  systems  into  which  to 
introduce  the  nonlinearities.  These  were  chosen  on  the  bases  (l)  that 
the  frequency  at  which  the  direct  polar  plot  of  the  system  forward 
transfer  function  crosses  the  negative  real  axis  be  within  the  limits 
of  the  analog  computer  and  its  recording  devices,  so  that  the  predic¬ 
tions  can  be  checked  experimentally  without  a  lot  of  time  scaling, 
and  (2)  that  the  system  open  loop  transfer  function  be  of  higher  than 
second  order  so  that  the  describing  function  method  is  preferable  over 
the  phase  plane  method  of  analysis.^  (Also,  the  system  open  loop 
transfer  functions  must  be  of  higher  than  second  order  before  self- 


2 

A  function  is  said  to  be  odd  periodic  if  f(t}=— f(-t). 

A  function  is  said  to  be  odd  harmonic  if  f ( t) =— f ( t4T/2) ,  where  T  is 
the  period. 


^For  a  description  of  phase  plane  analysis  see  the  references  (Ref 
12:65ff). 


5 


GGC/EE/64-16 


sustained  oscillations  can  occur  with  nonlinearities  that  do  not 
introduce  a  phase  shift.)  Nine  different  linear  forward  transfer 
functions  to  be  used  in  'unity  feedback  systems  were  chosen  on  these 
bases.  They  are  described  qualitatively  in  Table  I  on  the  following 
page.  The  table  also  indicates  the  nonlinearities  with  which  each 
was  coupled.  Since  there  are  to  be  five  describing  functions  inves¬ 
tigated  for  each  nonlinearity,  the  scope  of  the  problem  involves  220 
amplitude  and  frequency  predictions,  and,  since  each  system  must  be 
checked  on  the  analog  computer  as  a  basis  for  evaluation  of  the  accu¬ 
racy  of  the  predictions,  44  analog  computer  runs. 

Analysis  of  the  Problem 

Assumptions.  The  basic  assumptions  which  one  must  make  in  order 
to  use  the  describing  function  method  of  analysis  to  determine  the 
amplitude  and  frequency  of  self-sustained  oscillations  are  (l)  that 
the  system  is  unforced  and  does  not  vary  with  time,  (2)  that  the  non¬ 
linearity  is  separable  and  is  not  frequency  dependent,^  and  (3)  that 
the  linear  transfer  function  provides  sufficient  low-pass  filtering 
to  permit  excluding  the  higher  harmonics  from  consideration  (Ref  2:348). 
These  assumptions  are  quite  good  for  the  nonlinearities  and  linear  sys¬ 
tems  to  be  studied,  and  they  are  also  generally  true  in  practice. 

Since  these  assumptions  are  well  understood  by  users  of  the  describ¬ 
ing  function  method,  they  will  not  be  further  belabored. 


4 

Some  authors  extend  the  describing  function  to  frequency  dependent 
nonlinearities  as  well;  however,  as  indicated  in  the  definition  on 
page  1,  it  would  then  be  a  "generalized  transfer  function"  not  a 
"describing  function." 
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Table  I 

Qualitative  Description  of  the  Nine  Linear 
Forward.  Transfer  Functions  to  Be  Studied. 


Linear  Forward 

Transfer  Functions 

Nonlinearities 

Applied 

Type* 

Order** 

Stability 

Appro xi sate 
(Direct  Polar 
Plot)  Heal  Axis 
Crossing*** 

0 

3 

Unstable 

-  3\ 

1 

3 

Unstable 

-  3  \ 

2 

3 

Unstable 

-  3  1 

Saturation,  Satura- 

tion  with  Dead  Zone, 

Ideal  Relay,  Relay 

2 

3 

Unstable 

-10  > 

with  Hysteresis, 

1 

Relay  with  Dead  Zone, 

0 

5 

Unstable 

-10  I 

Relay  with  Dead  Zone 

and  Hysteresis. 

1 

5 

Unstable 

-10  / 

2 

5 

Unstable 

-10  / 

Conditionally 

r\ 

Stable  # 

~  3  \ 

> 

Dead  Zone 

\ 

\ 

Conditionally 

10  / 

Stable  # 

-loy 

♦Type  0,  1,  and  2  systems  are  wiled  position,  velocity,  and 
acceleration  control  systems,  respectively,  by  some  authors.  The 
numbers  0,  1,  and  2  refer  to  the  pover  of  the  separable  s,  the 
Laplacian  operator,  in  the  denominator  of  the  linear  transfer 
function. 

♦♦The  number  of  poles  minus  the  number  of  zeros. 

♦♦♦The  gain  of  each  system  was  set  for  one  of  the  two  distinct 
crossings  listed  in  order  to  provide  a  vide  variation  in  the  ampli¬ 
tudes  of  the  sustained  oscillations. 

#  It  is  necessary  to  use  a  conditionally  stable  system  in  order  to 
reach  a  stable  limit  cycle  with  a  dead  zone  nonlinearity. 
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The  new  rms  describing  function  has  not  been  extended,  as  yet, 
to  the  analysis  of  nonlinearities  with  memory  (Ref  2:384)-  In  the 
author's  opinion,  however,  this  scheme  can  be  applied  quite  easily 
to  nonlinearities  with  memory  that  result  in  a  phase  shift  of  the 
output  wave  that  can  be  calculated  algebraically,  such  as  for  the 
hysteresis  nonlinearities  included  in  this  paper.  If  the  axes  are 
shifted  through  the  proper  phase,  the  output  of  the  nonlinearity 
becomes  an  odd-harmonic,  odd-periodic  function,  and  the  new  rms 
describing  function  method  can  easily  be  applied  to  determine  the 
magnitude  of  the  describing  function.  Furthermore,  the  shift  in 
phase  becomes  the  phase  angle  of  the  resulting  new  rms  describing 
function. 

The  idea  behind  the  corrected-conventional  describing  func¬ 
tion  is  that  perhaps  the  effect  of  harmonics  higher  than  the  first 
in  the  Fourier  Series  expansion  of  the  output  of  the  nonlinearity, 
each  receiving  a  different  degree  of  attenuation  through  the  linear 

portion  of  the  system,  can  be  approximated  by  including  the  unat- 

5 

tenuated  effect  of  the  third  harmonic  only.  This  effect  is  "aver¬ 
aged  in"  by  the  same  method  one  would  use  to  obtain  the  combined 
elfective  rms  value  of  two  sine  waves  of  different  frequencies  when 
the  rms  value  of  each  is  known.  The  square  root  of  the  sum  of  the 
squares  of  the  first  and  third  Fourier  harmonics  is  therefore  taken 
as  the  amplitude  of  the  equivalent  sinusoidal  output  of  the  non¬ 
linearity,  rather  than  just  the  first  harmonic  as  in  the  cor.ven- 


^This  explanation  assumes  that  the  output  of  the  nonlinearity  is 
odd  harmonic  in  which  case  the  third  harmonic  is  the  first  har¬ 
monic  above  the  fundamental. 
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tional  describing  function.  It  would  seem  that  the  corrected- 
conventional  describing  function  would  tend  to  overcompensate  for 
the  error  caused  by  higher  harmonics,  especially  as  the  order  of  the 
system  is  increased.  However,  the  magnitude  of  the  corrected-conven- 
tional  describing  function  is  generally  less  than  that  of  the  new 
rms  describing  function,  thus  making  it  applicable  to  a  wider  range 
of  intermediate-order  systems  than  the  new  rm3  describing  function 
which  tends  to  overcompensate  (Ref  3 *1321). 

Criteria.  It  would  be  nice  if  the  engineer  could  be  guaren- 
teed  that  one  of  the  describing  function  schemes  could  provide  him 
with  predictions  that  were  always  within  ten  percent  of  the  actual 
value  and  generally  within  three  percent  of  the  actual  value  of  the 
amplitude  and  frequency  of  the  self-sustained  oscillations.  To 
prove  rigorously  that  a  certain  describing  function  generating 
scheme  would  always  provide  describing  functions  which  were  within 
ten  percent  of  the  true  outcome  would  certainly  be  a  complex  and 
difficult,  if  not  impossible,  task.  It  is  possible  to  postulate, 
however,  based  on  the  more  or  less  arbitrary  selection  of  linear  sys¬ 
tems  for  this  experiment,  that  a  particular  method  generally  pro¬ 
vides  predictions  which  are  within  three  percent  of  the  true  outcome, 
and/or  is  generally  better  than  the  other  methods  of  prediction. 

Procedure.  The  preceding  criteria  intrinsically  call  for  a 
comparison  of  some  kind  between  the  various  describing-function 
generating  schemes.  This  was  done  in  the  following  manner:  (l) 
each  of  the  35  describing  functions  previously  mentioned  was  derived, 
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(2)  the  direct  polar  plots  of  the  linear  systems  and  the  negative 
inverse  plots  of  the  describing  functions  were  compared  to  determine 
the  intersection  points,  (3)  the  stable  intersections  were  inter¬ 
preted  as  predictions  of  the  amplitude  and  frequency  of  self-sus¬ 
tained  oscillations  in  the  unforced,  closed-loop  system  employing 
the  nonlinearity,  (4)  the  predictions  were  tabulated  along  with 
experimental  data  obtained  by  simulation  of  the  actual  nonlinear 
system  on  the  analog  computer,  and  (5)  the  experimental  results  were 
used  as  a  basis  of  comparison  of  the  accuracy  of  the  various  describ¬ 
ing  functions  so  that  a  preferable  describing-function  generating 
scheme  could  be  selected. 

By  this  procedure  a  "best"  method  of  describing  each  of  the 
selected  nonlinearities  will  be  postulated.  In  this  way  some  guide¬ 
lines  will  be  added  to  the  describing  function  method  of  nonlinear 
circuit  analysis  so  that  the  control  engineer  may  obtain  quanti¬ 
tative  as  well  as  qualitative  results  from  his  work.  In  addition, 
the  various  methods  of  transient-response  prediction  which  depend 
on  the  accuracy  of  the  describing  function  may  be  improved.^ 


For  example,  computing  the  closed-loop  frequency  response  of  non¬ 
linear  systems  as  described  in  the  references  (Ref  12:195-208). 
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II.  Linear  Sya tarns 


Table  I  in  Chapter  I  qualitatively  describes  the  forward  trans¬ 
fer  functions  of  the  linear  systems  selected  as  representative  of 
those  often  encountered  in  practice  for  the  purpose  of  study.  Since 
the  conclusions  and  recommendations  postulated  in  the  last  chapter  of 
this  paper  are  largely  dependent  on  just  how  arbitrary  and  represent¬ 
ative  the  linear  systems  are,  this  entire  chapter  is  devoted  to  their 
description. 


General  Block  Diagrams 

The  general  block  diagram  of  the  linear  systems  discussed  in 
the  introduction  is  shown  in  Figure  1(a)  below.  Figure  1(b)  shows 
the  same  system  with  a  nonlinearity  introduced.  The  nomenclature 
indicated  in  Figure  1(b)  will  be  used  throughout  the  remainder  of 
this  paper. 


wmwmm 

Mam  B  j  Ui  —Ha 

\  — 

* 

r_ 

(a) 


Figure  1 


oo 


Unity  feedback  systems:  (a)  linear,  (b)  nonlinear,  where  G  ■  linear 
forward  transfer  function,  S  «  describing  function  which  represents 
the  nonlinearity  when  x(t)  -  X  sin  cot. 
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Linear  Transfer  Functions 

Three  third-order  and  three  fifth-order  linear  forward  transfer 
functions  were  originally  chosen  to  he  used  in  this  study.  Their 
transfer  functions  are  as  follows,  each  written  in  two  common  forms: 


180 


°1(8)  -  (  8+l)  (  8+2}  (s+3)"  or 

°1^  “  (l+j^)  ( 1+jO^  50)  ( 1+ jO.  333^)  ^ 


G2(«) 


18 

s(s+l) (s+2) 


or 


g2(o^) 


1  +3W)  ( 1+jO.  5*o) 


(3) 


03(b) 


4.76(b+0.1) 

s2(s+l)2 


or 


G3(j«) 


0.476(l+nl0«o) 

(at«)2(i+o<^)2 


(4) 


°4^S^  *  ( s+l )( s+2 )( S+3 )( s+4 )( 8+57 


or 


°4^  "  T 1+H ( 1+jO.  y>)  ( 1+jO.  333w)  ( 1+jO.  2»  C 1+ jO. 2“') 


r  301 _ 

y  m  s( s+l) ( s+2) (s+3) ( S+4)  or 

f  \  1 2  54^ 

G^Co^)  ju(l+jco)'(i+j0.5«o)'(  i+j0.333‘*’)  ( 1  + j0.2^*j) 


(6) 


g6(s)  ,  .-^i^esto_,i  1 

s2(s+l)2(s+5)2 


or 
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0.97i7(  1+jlOto) _ 

(  ja?)2(l+jco)2(l+j0.2  to)2 


In  a  block  diagram  arrangement  such  as  that  of  Figure  1(a),  E’qs 
(2),  (3),  and  (4)  become  type  G,  1,  and  2  systems,  respectively. 

They  are  third  order  because  their  denominators  are  of  order  three 
higher  than  their  numerators.  Furthermore,  their  gains  have  been 
adjusted  so  that  they  cross  the  real  axis  of  their  direct  polar  plots 
at  approximately  -3.  (See  Figures  2,  3,  and  4  at  the  end  of  this 
chapter.)  Similarly,  Eqs  (5),  (6),  and  (7)  become  type  0,  1,  and  2 
systems,  respectively.  They  are  fifth  order  because  their  denomi¬ 
nators  are  of  order  five  higher  than  their  numerators.  Their  gains 
have  been  adjusted  so  that  they  cross  the  real  axis  of  their  direct 
polar  plots  at  approximately  -10.  (See  Figures  5,  6,  and  7.  at  the 
end  of  this  chapter.)  These  transfer  functions  are  believed  to  be 
somewhat  representative,  since  they,  as  most  well-designed  control 
systems,  have  no  more  than  three  to  five  dominant  (closed-loop) 
roots,  and  also,  the  three  types  of  systems  generally  encountered 
in  practice  are  included.  Furthermore ,  by  establishing  two  sep¬ 
arate  real  axis  crossing  points  for  each  order  group,  both  "large" 
and  "small"  (relatively  speaking)  limit  cycle  amplitudes  can  be 
obtained  when  the  nonlinearities  are  introduced.  However,  in  order 
to  determine  if  the  difference  in  accuracy  between  the  predictions 
ior  G^(s),  G^s),  and  G^(s)  ajid  the  predictions  fcr  G^(s),  G^(s), 
and  G.i s)  was  due  to  limit  cycle  amplitude  or  system  order,  the 
gain  ol  G^(s)  in  Eq  (4)  vas  increased,  allowing  it  to  cross  the  real 
axis  of  the  direct  polar  plot  at  the  same  point  as  Gg(s)  in  Eq  (7). 
Thus  a  seventh  linear  transfer  function  was  introduced  for  study: 
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(») 


16.57(840.1) 

.2(«f  D2 

V3W) 


or 

,  1.657(1+310  «o) 

(j^)2(l4-i  w)2 


(8) 


The  direct  polar  plot  of  this  function  is  similar  to  that  of  Figure 
4,  except  that  the  real  axis  crossing  occurs  at  -10.339  instead  of 
-2.975- 

Eqs  (2)  through  (8)  comprise  a  set  of  representative  linear 
functions  into  which  to  introduce  the  nonlinearities.  Upon  analysis 
of  these  systems ,  however,  it  was  found  that  none  vould  produce 
stable  oscillations  with  a  dead  zone  nonlinearity.  Therefore,  in 
order  to  study  the  accuracy  with  which  each  type  of  dead  zone  de¬ 
scribing  function  can  predict  the  amplitude  and  frequency  of  self- 
sustained  oscillations,  it  is  necessary  to  first  devise  a  linear 
system  which  will  produce  stable  oscillations  when  dead  zone  is 
introduced.  Eqs  (9)  and  (10)  below  are  the  forward  transfer  func¬ 
tions  of  two  such  systems. 


0.(0) .  afegsteal. 

7  s3(s+10)2 

0?(  j  “>) 


or 

2.46626(1+3  v)2 

(jw)' 3  (1+30.  l^)2 


07A(b) 


s3(s+10)2 


or 


g7^w) 


l  8.26756(1+;)  ^)2 

(3^)3(l+j0.1  co)2 


(9) 


(10) 


14 


ggc/ee/ 64-16 


As  can  be  seen  fron  an  analysis  of  Figures  8(a)  and  (b)  at  the 
end  of  this  chapter,  these  conditionally  stable,  type  3,  third-order 
linear  systems  will  exhibit  stable  oscillations  when  coupled  with  a 
dead  zone  nonlinearity.  (The  direct  polar  plot  of  G?.(s)  is  similar 
to  that  in  Figure  3(a) ,  except  that  the  real  axis  crossings  occur  at 
-0. 6o6o33  and  -10,  instead  of  at  -0.20605  and  -3  as  shown  for  system 
G7(s).) 

Eqs  (2)  through  (10)  are  the  forward  transfer  functions  of  the 
so-called  arbitrary,  representative  linear  systems  into  which  the 
nonlinearities  covered  in  this  study  were  introduced,  with  the 
exceptions  noted.  In  future  discussions  in  this  paper,  the  linear 


functions  of  Eqs  (2)  through  (10)  will  often  be  referred  to  by  their 


subscripted  symbols;  i.e.,  G^s)  or  G^A(ju>)),  as  has  been  done  in 
several  places  in  this  chapter.  With  the  linear  systems  thus  estab¬ 


lished  and  understood,  the  next  step  is  to  generate  some  describing 


functions  for  the  nonlinearities  mentioned  in  Chapter  I.  This  will 


be  done  in  the  next  chapter. 
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polar  plot  of  Sq 
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Figure  6 


Direct  polar  plot  of  Eq  (6) 
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III.  Describing  Functions 

Tbs  general  definition  of  a  describing  function  is  given  on 
page  one  in  the  introduction  of  this  paper.  In  this  chapter  that 
definition  is  expanded  in  order  to  define  each  of  the  specific 
describing-function  generating  schemes  included  in  the  study:  con¬ 
ventional,  minimum  average  error,  equivalent  gain,  new  rms,  and 
corrected-conventional.  Then  the  describing  function  that  each 
scheme  produces  is  listed  for  each  of  the  nonlinearities  included 
in  the  study:  dead  zone,  saturation,  saturation  vith  dead  zone, 
ideal  relay,  relay  vith  hysteresis,  relay  vith  dead  zone,  and  relay 
with  both  dead  zone  and  hysteresis. 

Definitions 

Since  the  describing-function  generating  schemes  differ  only  in 
the  vay  the  equivalent  output  sine  vave  from  the  nonlinearity  is 
chosen,  the  definition  of  each  scheme  can  be  placsd  in  a  more  or  less 
standardized  format  based  on  the  general  definition  of  the  describ¬ 
ing  function.  This  procedure  vill  facilitate  comparison  of  the 
methods.  In  each  of  the  following  definitions 

x(t)  *  X  sin  *ot  =  the  input  to  the  nonlinearity} 

y(t)  =  the  actual  nonsinusoidal  output  from  the  nonlinearity. 

The  Conventional  Describing  Function.  The  conventional  describ¬ 
ing  function  is  equal  to  the  amplitude  and  phase  of  an  equivalent 
output  sine  vave,  chosen  as  the  fundamental  frequency  term  in  the 
Fourier  series  expansion  of  the  actual  nonsinusoidal  output  of  the 
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nonlinearity,  divided  by  the  amplitude  and  phaee  of  the  input  eine 
wave.  In  equation  fora  (Bef  12*146) 


(A^+JL2)*  /tan”1  A  * 

i- _ 


X  /0 


where 


3  r. 

TT  >0* 


y(  t)  coe  o>t  d  a>t 


B1  “  8in  d 


(11) 


The  Minimum  Average  Error  Describing  Function.  The  minimum 
average  error  describing  function  is  equal  to  the  amplitude  and  phase 
of  an  equivalent  output  sine  wave,  so  chosen  as  to  minimize  the  aver¬ 
age  error  between  itself  and  the' actual  nonsinusoidal  output  of  the 
nonlinearity,  divided  by  the  amplitude  and  phase  of  the  input  sine 
wave  (Ref  2:361).  Let  be  the  phase  difference  between  the  input 
sine  wave  and  the  output  periodic  waveform.  (The  phase  of  the  out¬ 
put  is  taken  as  the  shift  in  axes  necessary  to  make  the  output  wave  odd 
periodic.)  The  average  error  can  be  minimized  by  setting  the  dif¬ 
ference  of  the  area  under  a  half  cycle  of  the  actual  output  and  a 
half  cycle  of  the  equivalent  sinusoidal  output  equal  to  zero.  In 
equation  fora: 


average  error 


CTT+4>  (it+4> 

0  »  y(t)  d  an  -j ^  sin(  cot  -  <f>)  dwt 

(Tf+4> 

J<p  dwt  “  2Yi  ( 12) 

where  8in(tut  -/)  is  the  equivalent  output  sinusoid.  From  Eqs 
v 12)  and  (l)  the  describing  function  is,  ther tf'avn 5 

1  Cl T+<j> 


N  = 


y(t)  dcot 


(13) 
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The  hquivalent  Gain  Describing  Function.  The  equivalent  gain 
describing  function  is  equal  to  the  amplitude  and  phase  of  an  equiv¬ 
alent  output  sine  wave,  chosen  to  have  an  amplitude  equal  to  the  pealc 
value  of  the  actual  nonsinusoidal  output  of  the  nonlinearity,  divided 
by  the  amplitude  and  phase  of  the  input  sine  wave  (Ref  10:217).  In 
equation  form 


where  ^  is  the  phase  difference  between  the  input  sine  wave  and  the 
output  periodic  waveform.  (The  phase  of  the  output  is  taken  as  the 
shift  in  axes  necessary  to  make  the  output  wave  odd  periodic.) 

The  New  RMS  Describing  Function.  The  new  rms  describing  function 
is  equal  to  the  amplitude  and  phase  of  an  equivalent  output  sine 
wave,  so  chosen  as  to  have  the  same  rms  value  as  the  rms  value  of  the 
actual  nonsinusoidal  output  of  the  nonlinearity,  divided  by  the 
amplitude  and  phase  of  the  input  sine  wave.  In  equation  form  (Ref 
2:382) 


where  f  is  the  phase  difference  between  the  input  sine  wave  and  the 
output  periodic  waveform.  (The  phase  of  the  output  is  taken  as  the 
'••Vi ft  in  fives  necessary  to  make  the  output  wave  odd  periodic.) 
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The  Corrected-Conventional  Describing  Function.  The  corrected- 

conventional  describing  function  is  equal  to  the  amplitude  and  phase 

of  an  equivalent  output  sine  wave,  chosen  to  have  an  amplitude  equal 

to  the  square  root  of  the  sum  of  the  squares  of  the  coefficients  of 

the  fundamental  and  third  harmonic  in  the  Fourier  series  expansion  of 

7 

the  actual  nonsinusoidal  output  of  the  nonlinearity'  and  a  phase 
equal  to  the  phase  of  the  fundamental  frequency  term  in  the  expansion., 
divided  by  the  amplitude  and  phase  of  the  input  sine  wave.  In 
equation  form 


5 


where 


(16) 


2ftr 

Jo 


®1  "tt\  y(t)  sin  tut  d  cot 


2 

#  -  *  rf!  l(t)  008  3  Wt  d  cot 

J  '■  Jo 


B 


3 


sin  3  cot 


dcot 


7 

'Only  nonlinearxties  which  produce  odd-harmonic  outputs  are  considered 
thus  far  in  the  definition*  therefore,  the  third  harmonic  is  the  first 
harmonic  above  the  fundamental. 
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The  application  of  the  5  foregoing  definitions  to  the  7  non- 
linearities  to  be  considered  yields  the  35  describing  functions  to 
be  studied. 

Derivations 

The  35  describing  functions  are  listed  in  Tables  II  through  VIII 
on  the  following  pages,  with  one  table  for  each  nonlinearity  to  be 
considered.  The  various  describing  functions  listed  are  nondimen- 
sionalized  to  the  extent  necessary  to  enable  graphic  presentation  of 
describing  function  by  a  single  curve  or  family  of  curves.  These 
describing  functions  were  derived  by  applying  Eqs  (ll),  (13),  (14), 
(15),  and  (16)  to  each  nonlinearity.  A  list  of  where  some  of  these 
derivations  can  be  found  is  given  in  Appendix  E.  The  derivations  of 
the  corrected-conventional  describing  functions,  which  are  new,  can 
be  found  in  Appendix  A,  and  plots  of  these  can  be  found  in  Appendix 
D.  Appendix  B  contains  the  derivations  of  the  new  rms  describing 
functions  for  a  relay  with  hysteresis  and  a  relay  with  both  dead  zone 
and  hysteresis.  These  describing  functions  had  not  heretofore  been 
derived  since  Gibson  had  not  yet  extended  the  new  rms  describing 
function  to  nonlinearities  which  are  not  single-valued  (Ref  2:384). 

With  the  linear  systems  of  Chapter  II  and  the  describing  func¬ 
tions  of  this  chapter  thus  established,  it  is  now  possible  to  proceed 
to  the  amplitude  and  frequency  predictions  and  the  comparison  of 
those  predictions  with  experimental  results. 
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Table  II 


Deaoribing  Functions  for  Dead  Zona  aa  Derived  by  Varioue  Deeoribing- 

Punotion  Generating  Sohemes. 

For  each  describing  funotion  In  the  tablet 

Definition  of  nonlinearity  and  paraaeterst 

Definition  of  variablet  -  ain“*  -gjr 

Conditional*  "^2  Inputt 


Deaoribing-Funotion 
Generating  Scheme 

Vondimenaionalized  Deaoribing  Funotion 

Conventional 

|  -  frr-  2*±  -  .in  2^) 

Minimum  Average 

Srror 

1  -  <**  -i),ia  ad +  °°,ocd 

Equivalent  Gain 

|  -  1  -  am  «a 

In  HKS 

Correoted- 

Gonventional 

f  .  i  j0T-  2«d  -  .in  2«4)2  +  (-|  .in  4*d  * 

a  r  —.A 

■|  ain  2^  sln2o^  -  ^  ain 

♦This  condition  implies  that  d/2*£  X<Coo.  Also,  if  X^d/2, 
there  is  no  output  and  N/K  -  0. 
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Table  III 

Describing  Functions  for  Saturation  as  Derived  by  Various  Describing- 

Function  Generating  Sohemes. 

For  each  desoribing  function  in  the  tablet 

Definition  of  nonlinearity  and  parameters 

—1  s 

Definition  of  Variable  t  OLg  m  Bin 

Conditions^  0^  ^^.TC/2  Inputt  z  -  X  sin  tut 


put put 


Desoribing-Funotion 
Generating  Scheme 

Mondimensionalised  Desoribing  Function 

Conventional 

Minimum  Average 

Error 

^  -  1  +  sin.  <Xm  -  cos 

Equivalent  Gain 

|  "  8in  «• 

ffev  RMS 

i .  [? . 

Correoted- 

Conventional 

^  ^  |(2«^  +  8in  2wm)2  +  ("3  Bin-  2l%  ~ 

■|  sin  40^  -  sin  2 sin2o<^) 

♦This  condition  implies  that  S/2^.X<oo.  Also,  if  X  ^  S/2, 
there  is  a  linear  output  and  n/k  «  1. 
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Table  IV 


Describing  Funotion*  for  Saturation  with  Dead  Zone  as  Derived  through 
the  Application  of  Various  Describing-Function  Generating  Scheme*. 


For  eaoh  describing  funotion  in  the  tablet 

Definition  of  nonlinearity  and  parameters t 


t  output 


S 

2 


d 

'Z 


Definition  of  variables t  (X,  -  sin 


-1  d 


21 

sin”1^ 


d 

2 


Conditional*  Od  OL.d.  OL  4=  77/2  Input t  x  -  X  simcDt 

CL  • 


s 

2 


Desoribing-Funotion. 
Generating  Sohene 

Hondimensionalised  Describing  Funotion 

Conventional 

|  "  ?<2*,  -  2^  +  sin  2c^  -  sin  2^) 

Minimum  Average 

Error 

e»  wr* 

^  -  (*4  -  oc^sin  <x.d  +  (tj  -  o^)(8in.  —  sin  et4)  + 

oos  «d  -  oos 

Equivalent  Gain 

~  -  sin  (Xm  -  sin  <*4 

Mew  HMS 

1  "  (?{2<*a  '  *d>  -  2^  -  3  sin  2^  + 

S  jobs  +  («B  -  sin  p^jsine<d  + 

4(f  -  «d)  «ia-2K4  +  4(f  -  «,) 

Correoted- 

Conventional 

1  "  ?£2*b  "  2*d  +  8in  2<xs  “  8in  2<^)2  + 

(■^  sin  2*<b  -  ^  sin  2*4  -  sin  4<*B  + 

•|  sin  4<*d  -  *y  sin?c<B  ooe  + 

•y  sinV^  oos 

»These  conditions  imply  that  S/2  4X<oo  and  S  >■  4.  If  X  <  S/2, 
Table  II  applies. 
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Describing  Functions  for  an  Ideal  Relay  aa  Derived  through  the 
Application  of  Various  Desoribing-Funotion  Generating  Scheme*. 


For  each  describing  function  in  the  table j 

Definition  of  nonlinearity  and  parameters:  - 

Conditions*  OZ  X  Zl  oo  Input:  x  -  X  sin  u)X 


loutput 


X 


-X 


input 


Describing-Function 
Generating  Scheme 

Nondimensionalised  Describing  Function 

Conventional 

* .  ± .  i&gia 
x  7m  x 

Minimum  Average 

Error 

N  7T  1. 571 

K  "  ZX  “  X 

Equivalent  Gain 

N  1 

K  "  X 

Kev  RMS 

K  VT  1.414 

K  "  X  "  X 

Corrected- 

Conventional 

N  4  ./10  1.342 

x  "imi  9  “  x 

3] 
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Table  VI 

Describing  Functions  for  a  Relay  vith  Hysteresis  as  Derived  through 
the  Application  of  Various  Desorlbing-Funotlon  Generating  Soheoes. 

For  eaoh  describing  function  in  the  tablet 


Definition  of  nonlinearity  and  parameters! 


Definition  of  variable t 


-1  h 


*h-Bin  -5 


output 

_ 

h 

2 

h 

,  2 

input 

_  f 

Conditions*  h/2 4^  XZ!oo  Input*  x  -  X  sin  u>t 


Desoribing-Funotion 
Generating  Scheme 


Conventional 


Minimum  Average 
Error 


Equivalent  Gain 


Rev  RMS 


Correoted- 

Conventional 


Vondimensionalised  Describing  Function* 


» .  hm 

K  I 


3  .  1-571  .*3“n 
K  X 


I  1  "3Wh 
K  "  X 


3  .  i^i! 

K  X 


I  .  .H243 

K  X  * 


*  "j0ih 

e  is  complex  number  notation  for  1  /  n. 
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Table  VII 


Describing  Punotione  for  a  Relay  vith  Dead  Zone  ae  Deri-red  through 
the  Application  of  Various  Describing-Function  Generating  Schemes. 

For  eaoh  describing  funotion  in  the  table* 


Definition  of  nonlinearity  and  parameters* 

-1  d 


Definition  of  variable* 


“  sin 


2X 


Conditions**  0^.  7T/2  Input* 


Desoribing-Function 
Generating  Scheme 

fiondimensionalized  Describing  Funotion 

Conventional 

£4  -  !•  .In  2*d 

Minimum  Average 

Error 

If  -  (7T-  2<Kd)Blno^ 

Equivalent  Gain 

f  -  2  •*»  «d 

Rev  EMS 

~  -  2  ^2  -  -^)sin2(J^ 

Correoted- 

Conventional 

T  "  0in  2*a£  +  “  4  »i»2*d)*] 

♦This  condition  implies  that  d/2  ^  X<oo.  Also,  if  X  ^  d/2, 
there  is  no  output  and  Kd/^C  -  0. 
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Table  VIII 


Describing  Functions  for  a  Belay  vith  Dead  Zone  and  Hysteresis  as 
Derived  through  the  Application  of  Various  Describing-Function 

Generating  Schemes. 

toutjait 


For  eaoh  describing  funotion  in  the  table: 

Definition  of  nonlinearity  and  parameters: 
Definition  of  Variables:  iX  -  sin-1(4±r) 


T  f 

h — d- 

— > 

M 

1 

d _ .input 

I 

-ili 

/3  -  7T-  sin"1(~r) 


Conditions:*  oC  Input:  x  -  X  sinwt 


De s oribing-Fun ot i on 
Generating  Scheme 

Hondimensionalised  Describing  Function 

Conventional 

it~R~  +&  \ 

s&sL.£.in«  .  2'  2  5 

Kind mum  Average 

Error 

,ijr  «.  +/$  \ 

.  2‘  2  ' 

Equivalent  Gain 

,j<2'  2  5 

Hew  HMS 

r~  i  +  0  \ 

hiail-a2^-*).!*2*]*  ,J<2~  2  1 

Correoted- 

Conventional 

4  '»‘Z(S~K)  j  . l*  2  } 

♦This  condition  implies  that  (d+h)/2  ^  X<oo.  Also,  if  X.  ^  (d+h)/2 
there  is  no  output  and  H  (d+h) 

K  "  °* 
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IV.  Amplitude  and  Frequency  Predictions 

In  order  to  preserve  continuity  in  the  discussion  to  follow, 
only  one  nonlinearity  will  he  discussed  at  a  time.  That  is,  for 
each  nonlinearity  the  technique  used  to  predict  the  amplitude  and 
frequency  of  the  self-sustained  oscillations  will  he  presented, 
followed  by  a  comparison  with  the  experimental  results  and  a  dis¬ 
cussion  of  a  few  conclusions  which  may  be  drawn  from  the  comparison, 
before  moving  on  to  another  nonlinearity.  Thus  this  chapter  repre¬ 
sents  several  small  studies  within  a  larger  study,  but  this  is  a 
necessary  consequence  of  the  nature  of  the  problem.  The  answer  to 
the  question  "Which  describing  function  is  best?"  may  depend  on  the 
particular  nonlinearity  under  consideration. 

In  the  analysis  to  follow,  the  "gain,"  K,  of  the  nonlinearity 
is  sometimes  other  than  unity.  However,  this  does  not  defeat  the 
purpose  of  having  two  separate  axis  crossings  for  the  linear  portion 
of  the  system  (See  Figure  2  through  8  in  Chapter  II)  hy  changing 
the  loop  gain,  as  one  might  think.  The  main  reason  for  having  two 
separate  axis  crossings  for  the  linear  transfer  functions  was  to 
provide  an  adequate  separation  of  the  systems  into  two  groups,  so 
that  the  "gain,"  K,  of  each  nonlinearity  could  then  be  adjusted  to 
provide  a  "small"  amplitude  limit-cycle  group  and  a  "large"  ampli¬ 
tude  limit-cycle  group.  That  is,  exactly  where  the  "large"  signal 
and  "small"  signal  portions  of  the  describing  function  curve  occur 
depend  on  the  nature  of  each  individual  nonlinearity.  Therefore,  the 
value  of  K  must  be  selected  individually  for  each  nonlinearity,  so 
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that  the  separation  -will  actually  provide  what  nay  he  considered 
"large"  and  "small"  limit-cycle  amplitudes  for  that  nonlinearity. 

Saturation 

If  a  saturation  nonlinearity  is  introduced  into  one  of  the 
linear  systems,  G^s)  through  Gg(o),  or  G^(s),  in  the  manner  shown 
in  Figure  1(h)  on  page  11,  a  stability  analysis  will  reveal  that 
stable,  self-sustained  osoillations  will  result.  A  normalized  plot 
of  each  of  the  five  saturation  describing  functions  listed  in  Table 
III  on  page  29  is  shown,  in  Figure  9  on  pags  38.  The  values  of  the 
parameters  of  the  saturation  nonlinearity  were  ohosen  as  S  -  3  and 
K  -  0.5.  The  amplitude  and  frequency  predictions  by  each  of  the 
five  describing  functions  for  each  of  the  seven  systems  were  calcu¬ 
lated  in  the  following  manner: 

(1)  Use  Eqs  (2)  through  (8)  or  Figures  2  through  7  in  Chapter 

II  to  calculate  the  negative  real  axis  crossing  point  of  G  (jw)  on 

n 

the  direct  polar  plot.  The  value  of  Ga(jco)  at  the  negative  real 
axis  crossing  is  equal  to  -l/H,  and  the  value  of  cc?  at  the  negative 
real  axis  crossing  is  the  predicted  frequency  of  the  self-sustained 
oscillations. 

(2)  Calculate  the  value  of  -K/N  and  use  Figure  9  to  determine 
the  value  of  2X/S  for  each  system  and  each  describing  function. 

(3)  Calculate  the  value  of  X  from  the  known  values  of  2X/S  and 
S.  X  is  then  the  predicted  amplitude  of  the  self-sustained  oscil¬ 
lations  for  that  particular  describing  function  and  system. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Table  IX  on  page  39*  In  order  to  determine  the  accu- 
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racy  of  the  predictions,  each  of  the  seven  nonlinear  systems  was 
simulated  on  an  analog  computer.  (The  circuits  used  in  the  simula¬ 
tion  are  shown  in  Appendix  C.)  The  analog  computer  experimental 
results  are  also  shown  in  Table  IX,  along  with  the  percentage  of 
error  of  each  prediction  based  on  the  experimental  results. 

Table  X  on  page  40  is  a  comparison,  based  on  the  error  figures 
in  Table  IX,  of  the  various  saturation  describing  functions.  It 
shows  that  the  saturation  nonlinearity  seems  to  favor  the  conven¬ 
tional  describing  function,  although  it  can  be  seen  from  an  analysis 
of  Table  IX  that  the  accuracy  of  the  corrected-conventional  describ¬ 
ing  function  becomes  generally  better  as  the  amplitude  of  the  limit 
cycle  is  increased.  The  situation  is  somewhat  different,  however, 
when  dead  zone  is  added  to  the  saturation,  which  is  the  next  non¬ 
linearity  to  be  discussed. 
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Saturation  Combined  with  Dead.  Zone 

If  a  nonlinearity  consisting  of  saturation  combined  with  dead 
zone  is  introduced  into  one  of  the  linear  systems,  G^(s)  through 
G^(s),  or  G^a(s),  in  the  manner  shown  in  Figure  l(b)  on  page  11, 
a  stability  analysis  will  reveal  that  stable,  self-sustained  oscil¬ 
lations  can  exist.  A  normalized  plot  of  each  of  the  five  satura¬ 
tion  describing  functions  listed  in  Table  IV  on  page  30  is  shown 
in  Figure  10  on  page  43  for  the  case  of  S/d  =  2.  The  values  of  the 
parmeters  of  the  nonlinearity  were  chosen  as  S  =  4,  d  «  2,  and  K  *  1 
lor  the  purpose  of  analysis.  The  predictions  of  the  amplitude  and 
frequency  of  the  self-sustained  oscillations  by  each  of  the  five 
describing  functions  for  each  of  the  seven  systems  mentioned  were 
calculated  in  a  manner  similar  to  that  uBed  for  saturation: 

(1)  Use  Fqs  (2)  through  (8)  or  Figures  2  through  7  in  Chapter 
II  to  calculate  the  negative  real  axis  crossing  point  of  Gn(jo>)  on 
the  direct  polar  plot.  The  value  of  Gn( j w )  at  the  negative  real 
axis  crossing  is  equal  to  -l/tf,  and  the  value  of  U)  at  the  negative 
real  axis  crossing  is  the  predicted  frequency  of  the  self-sustained 
oscillations. 

(2)  Calculate  the  value  of  -K/K  and  use  Figure  10  to  determine 
the  value  of  2X/d  for  each  system  and  each  describing  function. 

(3)  Calculate  the  value  of  X  from  the  known  values  of  2X/d  and 
d.  X  is  then  the  predicted  amplitude  of  the  self-sustained  oscil¬ 
lations  lor  that  particular  describing  function  and  system. 

The  results  01  the  application  of  the  foregoing  procedure  are 
tabulated  in  Table  XI  on  page  44.  In  order  to  determine  the  accu¬ 
racy  ol  the  predictions,  each  of  the  seven  nonlinear  systems  was 
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simulated  on  an  analog  computer.  (The  circuits  used  for  the  simula¬ 
tion  are  shown  in  Appendix  C.)  'The  analog  computer  experimental 
results  are  also  shown  in  Table  XI,  along  with  the  percentage  of 
error  of  each  prediction  based  on  the  experimental  results. 

Table  XII  on  page  45  i8  a  comparison,  based  on  the  error  figures 
in  Table  XI,  of  the  various  saturation-plus-dead-zone  describing 
functions.  It  shows  that  the  saturation  +  dead  zone  nonlinearity 
Beems  to  favor  the  corrected-conventional  describing  function; 
furthermore,  analysis  of  Table  XI  shows  that  this  is  true  for  both 
large  and  small  limit-cycle  amplitude.  Thus,  the  addition  of  dead 
zone  to  the  saturation  must  have  increased  the  harmonic  content  of 
the  output  of  the  nonlinearity  to  the  point  that  a  correction  was 
needed  for  all  amplitudes  instead  of  just  for  high  amplitudes.  A 
similar  thing  happens  when  the  ideal  relay  is  analyzed. 
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Table  XI 

Prediction  of  the  Amplitude  and  Frequency  of  Self-Sustained  Oaoillationu  in  Seven  Different  Feedback 
Systems  with  Saturation  Combined  with  Dead  Zone  by  Varioue  Describing  Functions.** 
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Ideal  Relay 

If  an  ideal  relay  nonlinearity  is  introduced  into  one  of  the 
linear  systems,  Q^s)  through  G^(s),  or  G^(s),  in  the  manner  shovn 
in  Figure  l(h)  on  page  11,  a  stability  analysis  will  reveal  that 
stable,  self-sustained  oscillations  will  result.  A  plot  of  each  of 
the  five  ideal  relay  describing  functions  listed  in  Table  V  on  page 
31  is  shown  in  Figure  11  on  page  48.  The  value  of  K  was  chosen  to 
be  unity  lor  the  analysis.  The  predictions  of  the  amplitude  and 
frequency  of  the  self-sustained  oscillations  by  each  of  the  five 
describing  functions  for  each  of  the  seven  systems  mentioned  above 
were  calculated  in  the  following  manner: 

(1)  Use  Lqs  (2)  through  (8)  or  Figures  2  through  7  in  Chapter 

II  to  calculate  the  negative  real  axis  crossing  point  of  on 

the  direct  polar  plot.  The  value  of  Gn(j*>)  at  the  negative  real  axis 
crossing  is  equal  to  -l/lT,  and  the  value  of  at  the  negative  real 
axis  crossing  is  the  predicted  frequency  of  the  self-sustained 
oscillations. 

(2)  Calculate  the  value  of  -K/lT  and  use  Figure  11  or  the  describ¬ 
ing  function  equations  in  Table  V  to  determine  the  value  of  X,  the 
predicted  amplitude  of  the  self-sustained  oscillations  for  that 
particular  describing  function  and  system. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Table  XIII  on  page  49*  la  order  to  determine  the  accu¬ 
racy  of  the  predictions,  each  of  the  seven  nonlinear  systems  was 
simulated  on  an  analog  computer.  (The  circuits  used  for  the  simu¬ 
lation  are  shown  in  Appendix  C. )  The  analog  computer  experimental 
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results  are  also  shown  in  Table  XIII,  along  with  the  percentage  of 
error  of  each  prediction  based  on  the  experimental  results. 

Table  XTV  on  page  5°  i®  a  comparison,  based  on  the  error  figures 
in  Table  XIII,  of  the  various  ideal  relay  describing  functions.  It 
shows  that  the  ideal  relay  nonlineff.rity  strongly  favors  the  corrected- 
conventional  describing  function  for  the  systems  considered.  The 
same  thing  seems  to  be  true  of  all  the  various  relay  nonlinearities, 
insofar  as  amplitude  accuracy  is  concerned,  as  will  be  seen  in  the 
next  three  sections. 
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•All  of  the  describing  functions  listed  predict  the  same  frequency.  (Original  data) 

**K  =  1  (See  Table  V,  page  31). 
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Relay  with  Hysteresis 

If  a  nonlinearity  such  as  a  relay  with  hysteresis  is  introduced 
into  one  of  the  linear  systems,  G^(s)  through  Gg(s),  or  G^^(s),  in 
the  manner  shown  in  Figure  1(h)  on  page  11,  a  ataoility  analysis  will 
reveal  that  stable,  self-sustained  oscillations  will  occur.  A 
comparison  of  Tables  V  and  VI  (pages  31  and  32)  shows  that  the  magni¬ 
tude  variations  of  the  relay-with-hysteresis  describing  functions 
are  the  same  as  those  of  their  ideal  relay  counterparts;  however, 
the  hysteresis  introduces  an  amplitude-dependent  phase  shift.  Figure 


versus  I  for  the 


11  on  page  48  is  there! ore  a  valid  plot  of  -j 
relay  with  hysteresis.  In  order  to  analyze  the  relay-with-hystere¬ 
sis  describing  functions,  values  of  h  «  6  and  K  =  1  were  chosen, 
and  the  five  describing  functions  of  Table  VI  were  plotted  in  a 
psuedo  magnitude-angle  diagram  fashion  as  shown  in  Figure  12  on 
page  54-  (A  magnitude  rather  than  log-magnitude  scale  was  used  in 
order  to  preserve  accuracy.)  On  the  same  graph,  the  magnitude 
versus  single  relationship  of  each  of  the  seven  previously  mentioned 
linear  forward  transfer  functions  was  plotted.  The  intersection  of 
a  G  (jui)  curve  with  a  negative  inverse  describing  function  curve 
is  the  operating  point  at  which  that  particular  type  of  describing 
junction  predicts  a  stable  limit  cycle.  One  can  use  Figure  12, 
therefore,  to  make  the  amplitude  and  frequency  predictions  in  the 
following  manner: 


(l)  Use  the  magnitude -angle  condition  at  the  intersection  to 
solve  the  appropriate  equation  (Eqs  (2)  through  (8)  in 

Chapter  II)  lor<o.  This  value  of  w  is  the  predicted  frequency  cf 
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the  self-sustained  oscillations. 

(2)  The  magnitude  value  at  each  intersection  is  equal  to 
for  that  particular  system  and  describing  function.  Therefore 
Figure  11,  or  the  magnitude  portion  of  the  describing  function 
defining  equation,  can  be  used  to  determine  X,  the  predicted  ampli¬ 
tude  ol  the  self-su stained  oscillations. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Tables  XV  and  XVII  on  pages  55  and  57*  In  order  to 
determine  the  accuracy  of  the  predictions,  each  of  the  seven  non¬ 
linear  systems  was  simulated  on  an  analog  computer.  (The  circuits 
used  for  the  simulation  are  shown  in  Appendix  C.)  The  analog  com¬ 
puter  experimental  results  are  also  shovn  in  Tables  XV  and  XVII, 
along  with  the  percentage  of  error  of  each  prediction  based  on  the 
experimental  results. 

Table  XVI  on  page  $6  is  a  comparison  based  on  the  amplitude 
prediction  error  figures  in  Table  XV,  of  the  various  relay  +  hys¬ 
teresis  describing  functions.  It  shovs  that  the  relay-plus-hystere¬ 
sis  nonlinearity  seems  to  favor  the  corrected-conventional  describ¬ 
ing  function  insofar  as  amplitude  accuracy  is  concerned.  However, 
Table  XVIII  on  page  5 which  is  a  comparison  based  or.  the  frequency 
prediction  error  figures  in  Table  XVII,  shows  that  the  relay-vith- 
bysteresis  nonlinearity  favors  the  equivalent  gain  describing  func¬ 
tion  for  frequency  prediction  accuracy.  And,  from  an  analysis  of 
Table  XVII,  the  frequency  prediction  accuracy  of  the  equivalent 
gam  method  apparently  increases  as  the  limit-cycle  amplitude  is 
increased.  This  problem  of  a  variation  ir.  the  frequency  ore  fiction 
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Relay  with  Dead  Zone 

If  a  nonlinearity  such  as  a  relay  with  dead  zone  is  introduced 
into  one  of  the  linear  systems,  G^s)  through  G^(s),  or  G-^(s),  in 
the  manner  shown  in  Figure  1(b)  on  page  11,  a  stability  analysis 
will  reveal  that  stable,  self-sustained  oscillations  can  exist.  A 
normalized  plot  of  each  of  the  five  re lay-with-de ad-zone  describing 
functions  listed  in  Table  VII  on  page  33  ia  shown  in  Figure  13  on 
page  61.  For  the  analysis,  the  values  of  the  parameters  of  the  non¬ 
linearity  were  chosen  to  be  d  *  3  and  K  *  1.  The  predictions  of 
the  amplitude  and  frequency  of  the  self-sustained  oscillations  by 
each  of  the  five  describing  functions  for  each  of  the  seven  systems 
mentioned  were  calculated  as  follows: 

(1)  Use  Eqs  (2)  through  (8)  or  Figures  2  through  7  in  Chapter 
II  to  calculate  the  negative  real  axis  crossing  point  of  G  (jeo)  on 
the  direct  polar  plot.  The  value  of  Gq(o^)  at  the  negative  real 
axis  crossing  is  equal  to  -  l/K,  and  the  value  of  cj  at  the  negative 
real  axis  crossing  is  the  predicted  frequency  of  the  self-sustained 
oscillations. 

(2)  Calculate  the  value  of  -  K/Nd  and  use  Figure  13  to  deter¬ 
mine  the  value  of  2X/d  for  each  system  and  each  describing  function. 

(3)  Calculate  the  value  of  X  from  the  known  values  of  2X/d  and 

d.  X  is  then  the  predicted  amplitude  of  the  s elf-sustained  oscil¬ 
lations  for  that  particular  describing  function  and  system. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Table  XIX  on  page  62.  In  order  to  determine  the  accu¬ 
racy  of  the  predictions,  each  of  the  seven  nonlinear  systems  was 
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simulated,  on  an  analog  computer.  (The  circuits  used  lor  the  simu¬ 
lation  are  shown  in  Appendix  C.)  The  analog  computer  experimental 
results  are  also  shown  in  Table  XIX,  along  with  the  percentage  of 
error  of  each  prediction  based  on  the  experimental  results. 

Table  XX  on  page  63  i8  a  comparison,  based  on  the  error  figures 
in  Table  XIX,  of  the  various  relay-plus-dead-zone  describing  func¬ 
tions.  It  shows  that  the  relay  +  dead  zone  nonlinearity  seems  to 
favor  the  corrected-conventional  describing  function.  This  is  true 
because,  as  Table  XIX  shows,  for  small  amplitude  limit-cycles  when 
the  harmonic  content  'is  apparently  smaller,  the  conventional  and  the 
corrected-conventional  describing  functions  are  (nearly)  the  same; 
whereas,  when  the  limit-cycle  amplitude  is  increased,  the  corrected- 
conventional  describing  function  automatically  compensates  for  the 
increase  in  harmonic  content  of  the  output  of  the  nonlinearity. 

When  hysteresis  is  added  to  the  relay  with  dead  zone,  the 
problem  of  frequency  variation  between  the  predictions  again  arises. 
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Table  XIX 

Prediction  of  the  dmplitude  and  Frequency  of  Self-Sustained  Oscillations  in  Seven  Different  Feedback 
Systems  with  a  Relay  with  Dead  Zone  by  Various  Describing  Functions.** 
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of  the  Prediction  Errors  of  the  Various  Relay  +  Dead  Zone  Describing  Functions 


aCC/EE/64-16 


c 

o 

(0 

•H 

U 

a 

s 

o 

o 


<H  <S 


s 


d 


a  a 


►  s 


•  c 


an 

• 

« 

o. 

u 

hs 

o 

p 

Vt 

fi 

o 

o 

iti 

g 

<3 

*) 

r-4 

o 

CT*  P 

41 

c 

(H 

rH 

rH 

V 

cd 

E 

3 

3 

a * 

IQ 

V 

V 

(a 

rG 

0) 

p 

■H 

rH 

p 

On 

o 

(X 

'H 

cd 

( n 

£ 

•H 

O. 

03 

>> 

ri 

rH 

€> 

(d 

-P 

03 

§ 

•H 

rH 

r\ 

<n 

8 

c 

o 

® 

•H 

P 

{*> 

O 

o 

g 

G 

3 

V 

<H 

3 

o* 

tf) 

V 

c 

M 

•rl 

<H 

■ri 

t) 

f-i 

rG 

O 

P 

w 

4) 

D 

n 

f-l 

o 

e 

4h 

c 

9 

o 

H 

N 

e 

.G 

•g 

P 

4> 

ml 

§ 

+ 

>>  B 

flj 

« 

rH 

rH 

C 

£> 

U 

O 

%> 

cv 

-G 

P 

rH 

(4 

<w 

3 

o 

c 

*  -H 


63 


#  Indicates  the  minimum  magnitude  of  amplitude  error  in  that  row;  therefore,  the  column  with  the 
most#'s  is  the  column  of  what  is  poBsibly  the  best  describing  function  for  a  relay  with  dead 
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He lay  with  Dead  Zone  and  Hysteresis 

If  a  nonlinearity  such  as  a  relay  with  "both  dead  zone  and  hys¬ 
teresis  is  introduced  into  one  of  the  linear  systems,  G^(s)  through 
Gg(s),  or  G^a(s),  in  the  manner  shown  in  Figure  1(h)  on  page  11,  a 
stability  analysis  will  reveal  that  stable,  self-sustained  oscilla¬ 
tions  can  occur.  A  normalized  plot  of  the  magnitude  of  each  of  the 
five  relay-with-dead-zone-and-hysteresis  describing  functions  listed 
Table  VIII  on  page  34  is  shown  in  Figure  14  on  page  67  for  the  case 
of  d/h  -  1.  The  negative  magnitude  is  used  as  the  ordinate  to  permit 
comparison  with  the  other  describing  function  plots  in  this  chapter. 
The  angle  variation,  which  is  the  same  for  all  of  the  describing 
function  schemes,  is  shown  in  Figure  15  on  page  68.  In  order  to  ana¬ 
lyze  the  re  lay-with-dead-zone-and -hysteresis  describing  functions, 
values  of  d  ■  h  «  2  and  K  -  0.8  were  chosen,  and  the  five  describing 
functions  of  Table  VIII  were  plotted  in  another  psuedo  magnitude- 
angle  diagram  shown  in  Figure  16  on  page  69.  On  the  same  graph,  the 
magni tuie-versua-angle  relationship  of  each  of  the  seven  previously 
mentioned  linear  forward  transfer  functions  was  plotted.  The  inter¬ 
section  of  a  G^Cjto)  curve  with  a  negative  inverse  describing  func¬ 
tion  curve  is  the  operating  point  at  which  that  particular  type  of 
describing  function  predicts  a  stable  limit  cycle.  One  can  use  Fig¬ 
ure  16,  therefore,  to  make  the  amplitude  and  frequency  predictions  in 
the  following  manner: 

(l)  Use  the  magnitude -angle  condition  at  the  intersection  to 
solve  the  appropriate  Gn(jw)  equation  (Lqs  (2)  through  (8)  in  Chapter 
II)  for  w>.  Phis  value  of  «•->  is  the  predicted  frequency  of  the  self- 
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sustained  oscillations. 

(2)  The  magnitude  at  the  intersection  is  the  predicted  operating 
vali’.e  of  4  for  that  particular  system  and  describing  function. 

I  Xi  I 

Therefore,  Figure  14  can  be  used  to  determine  the  value  of  2X/d  by 
calculating  -|K/N(d+h)|  ana  then  picking  2X/d  off  of  the  proper  curve. 
Alternately,  the  angle  value  at  the  intersection  is  the  predicted 
operating  value  of  so  r'igure  15  can  be  used  to  find  2Z/d 

directly  from 

(3)  Calculate  the  value  of  X  from  the  knovn  values  of  2X/d  and 
d.  X  is  then  the  predicted  amplitude  of  the  self-sustained  oscilla¬ 
tions  for  that  particular  describing  function  and  system. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Tables  XXI  and  XXIII  on  pages  70  and  72,  respectively. 

In  order  to  determine  the  accuracy  of  these  predictions,  each  of  the 
seven  nonlinear  systems  was  simulated  on  an  analog  computer.  (The 
circuits  used  in  the  simulation  are  shown  in  Appendix  C.)  The  analog 
computer  experimental  results  are  also  shown  in  Tables  XX I  and  XXIII, 
along  with  the  percentage  of  error  of  each  prediction  based  on  the 
experimental  results. 

Table  XXII  on  page  71  is  a  comparison,  based  on  the  amputate 
prediction  error  figures  in  Table  XXI,  of  the  various  relay-with- 
dead-zone-and-hysteresis  describing  functionr .  It  shows  that  the 
relay  +  dead  zone  +  hysteresis  nonlinearity  favors  the  corrected- 
conventional  describing  function  insofar  as  amplitude  accuracy  is 
concerned.  However,  Table  XXIV.  on  yage  73,  which  is  a  comparison 
based  on  the  frequency  prediction  error  figures  in  Table  XXIII, 
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shows  that  the  relay  +  dead  zone  +  hysteresis  nonlinearity  favors 
the  equivalent  gain  describing  function  for  frequency  prediction 
accuracy.  Again,  as  in  the  case  of  a  relay  with  hysteresis,  the  accu¬ 
racy  of  the  equivalent  gain  frequency  prediction  seens  to  increase 
with  limit-cycle  amplitude  as  Table  XXIII  illustrates. 

The  last  nonlinearity  to  he  appraised  by  the  various  describing 
function  schemes  is  dead  zone.  Since  it  does  not  provide  a  stable 
limit  cycle  when  coupled  with  systems  G^(s)  through  Gg(s),  or  G^A(s), 
as  pointed  out  in  Chapter  II,  it  was  only  analyzed  with  two  condi¬ 
tionally  stable  systems,  or  rather,  the  same  system  at  two  gains. 
Therefore,  the  comparisons  are  less  conclusive  than  for  the  other 
nonlinearities,  but  they  are  still  interesting  and  necessary  for 
completeness. 
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Prediction  of  the  .Amplitude  of  Self -Sustained  Oscillations  in  Seven  Different  Feedback 
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#  Indicates  the  minimum  magnitude  of  error  in  that  row;  therefore,  the  column  with  the 
most#' a  is  the  column  of  what  is  possibly  the  best  describing  function  for  a  relay 
with  both  dead  zone  and.  hysteresis,  from  the  standpoint  of  amolitude  accuracy. 


Prediction  of  the  Frequency  of  Self -Sustained  Oscillations  in  Seven  Different  Feedback 
Systems  with  a  Relay  with  Dead  Zone  and  Hysteresis  by  Various  Describing  Functions.* 
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0.8  (bee  Table  VTII,  pa^e  34). 
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most  //  s  is  the  column  of  what  is  possibly  the  best  describing  function  for  a  relay 
with  both  dead  zone  and  hysteresis,  from  the  standpoint  of  frequency  accuracy. 
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Dead  Zone 

If  a  dead  zone  nonlinearity  ie  introduced  into  one  of  the 

linear  systems ,  G„(jto)  or  G_ .( j to) ,  in  the  manner  illustrated  in 

I  (A 

Figure  1(h)  on  page  11,  a  stability  analysis  will  reveal  that  sta¬ 
ble,  self-sustained  oscillations  will  result.  (See  Figure  8,  page 
22.)  A  normalized  plot  of  each  of  the  five  dead  zone  describing 
functions  listed  in  Table  II  on  page  28  is  shewn  in  Figure  17  on 
page  77*  Interestingly,  the  new  nns  describing  function  and  the 
corrected-conventional  describing  function  for  dead  zone  differ 
by  such  a  small  amount  that,  in  effect,  they  plot  as  the  same  curve. 
For  the  analysis,  the  amount  dead  zone  was  set  at  d  =  6  and  the 
gain  of  the  nonlinearity  at  K  =*  0.5.  The  predictions  of  the  ampli¬ 
tude  and  frequency  of  the  self-sustained  oscillations  by  each  of 
the  five  describing  functions  for  each  of  the  two  systems  men¬ 
tioned  were  calculated  as  follows: 

(l)  Use  Eqa  (9)  and  (10)  to  calculate  the  negative  real  axis 


crossing  point  of  the  direct  polar  plot  of  G7(jcu)  and  G..(  jw), 
respectively.  (Again,  see  Figure  8,  page  22.  The  pertinent  cross¬ 
ing  of  G^ 

discussed  in  Chapter  II.)  The  value  of  G^(ju))  at  the  crossing 
outside  of  the  -  1  point  is  equal  to  -  l/il,  and  the  value  of  co  at 
this  crossing  is  the  predicted  frequency  of  the  self-sustained 
oscillations. 

(2)  Calculate  the  value  of  -  k/n  and  use  Figure  17  to  deter¬ 
mine  the  value  of  2X/d  for  each  system  and  describing  function. 


(ju))  was  set  at  - 


and  G?A( 3  u ) 


course , 


(3)  Calculate  the  value  of  X  from  the  known  values  of  2X/d 
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and  d.  X  is  then  the  predicted  amplitude  of  the  self-sustained 
oscillations  for  that  particular  describing  function  and  system. 

The  results  of  the  application  of  the  foregoing  procedure  are 
tabulated  in  Table  XXV  on  page  76.  In  order  to  determine  the  accu¬ 
racy  oi  the  predictions,  both  of  the  nonlinear  systems  were  simu¬ 
lated  on  an  analog  computer.  (The  circuits  used  for  the  simula¬ 
tion  are  shown  in  Appendix  C.)  The  analog  computer  experimental 
results  are  also  shown  in  Table  XXV,  along  with  the  percentage  of 
error  of  each  prediction  based  on  the  experimental  results.  The 
frequency  error  for  system  G^^(ju))  was  far  greater  than  any  other 
frequency  error  in  the  entire  study.  This  can  probably  be  explained 
with  Figure  18  on  page  76.  As  can  be  seen  in  Figure  18,  the  basic 
describing  function  assumption  that  x(t)  is  approximately  sinus- 


it  was  a  good  assumption  in  every  other  case.  However,  it  is  inter¬ 
esting  to  note  that  the  amplitude  predictions  were  not  seriously 
affected  by  the  failure  of  this  assumption. 

Table  XXVI  on  page  79  is  a  comparison,  based  on  the  error 
figures  of  Table  XXV,  of  the  various  dead  zone  describing  functions. 
It  shows  that  the  dead  zone  nonlinearity  favors  both  the  new  rms  and 
the  corrected  conventional  describing  functions  since  they  give 
essentially  the  same  results  (to  three  or  four  place  accuracy). 

The  describing  function  comparison  tables  at  the  end  of  each 
section  of  this  chapter  have  used  four  error  criteria  as  a  basis 
for  drawing  conclusions  about  the  accuracy  of  each  type  of  describ¬ 
ing  function.  These  criteria  were  chosen  because,  un  the  author's 
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opinion,  the  "best”  describing  function  shouldt 

(1)  have  an  arithmetical  mean  error  (or  average  error)  very 
close  to  zero,  indicating  that  the  errors  encountered  with  that 
describing  function  vill  be  randomly  distributed  about  zero  error, 
rather  than  about  some  positive  or  negative  error  value,  so  that 
chances  of  having  an  error  of  almost  zero  are  maximum; 

(2)  have  a  median  error  close  to  the  value  of  the  arithmetical 
mean  error  to  indicate,  to  some  extent,  that  the  data  gathered  vas 
representative  and  not  "lopsided” 5 

(3)  have  the  smallest  maximum  magnitude  of  error  possible  in 
order  to  give  an  indication  of  its  maximum  "tolerance";  and 

(4)  have  the  smallest  arithmetical  mean  of  absolute  error  as 
an  indication  of  the  error  magnitude  that  can  usually  be  expected. 

This  concludes  the  analysis  of  the  describing  function  schemes 
as  they  apply  to  specific  nonlinearities.  In  the  next  chapter  an 
overall,  general  comparison  vill  be  made. 
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Figure  18 

■auataiaad  oa  dilation*  of  ay* taaa  (a)  0 7(ju)  aai  (t>)  O^Jw)  dth  daad 
_ _ _ _ _ alo— d  loop. 


Comparison  of  the  Prediction  Errors  of  tte  Various  Dead  Zone  Describing  Functions 
Prom  the  preceding  table  of  predictions  and  en*ors,  the  arithmetical  mean  error,  median  ■ 
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V.  Conclusions 


Some  readers  may  feel  that,  vhile  attempting  to  find  which 
describing  function  is  more  accurate,  the  author  has  only  compli¬ 
cated  the  problem  by  adding  still  another  scheme  to  the  field. 
However,  as  this  chapter  attempts  to  show,  although  the  corrected- 
conventional  describing  function  may  be  more  empirical  than  theo¬ 
retical,  it  does  remove  some  of  the  restraints  and  add  a  lot  of 
accuracy  to  the  analysis  of  a  wide  range  of  intermediate-order  non¬ 
linear  feedback  control  systems. 

Overall  Comparison  of  the  i)escribing-Function  Generating  Schemes 

■fable  XXVII  on  the  following  page  is  an  overall  comparison  of 
the  amplitude  prediction  errors  of  the  various  types  of  describing 
functions  included  in  this  study.  It  shows  quite  conclusively  that 
in  general,  to  obtain  amplitude  prediction  accuracy  when  using  the 
describing  function  technique  of  stability  analysis,  one  should 
use  the  corrected-conventional  describing-function  generating  scheme 
to  derive  the  describing  function.  Table  XXVII  shows  that  the  cor¬ 
rected-conventional  describing  function  is  generally  about  two 
percent  more  accurate  than  the  conventional  describing  iunction 
lor  intermediate-order  nonlinear  systems  such  as  those  included  in 
this  study.  It  also  shows  that,  in  general,  no  advantage  is  gained 
by  using  the  new  rms  describing  function  instead  of  the  conventional 
describing  function.  (However,  the  new  rms  describing  function  pro¬ 
vided  a  slight  improvement  lor  the  various  relay  nonlinearities 
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#  Indicates  the  minimum  magnitude  of  error  in  that  row;  therefore,  the  column  with  the 
most  r/' s  is  the  column  of  what  is  possibly  the  best  describing— function  generating 
scheme  to  use,  in  general,  in  order  to  attain  amplitude  accuracy. 
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discussed  in  Chapter  IV,  hut  the  improvement  vas  not  nearly  as 
significant  as  that  provided  hy  the  corrected-conventional  describ¬ 
ing  function.)  In  summary,  the  oorrected-oonventional  describing 
funotion  provides  certain  advantages  over  the  other  describing 
function  schemes  for  every  type  of  nonlinearity  tested  and,  in 
general,  it  has  a  mean  and  a  median  closer  to  zero  than  any  other, 
a  maximum  magnitude  of  error  less  than  any  other,  and  an  average 
absolute  error  less  than  any  other  describing  funotion  tested,  inso¬ 
far  a&  amplitude  prediction  accuracy  is  ooncemed. 

Trends  in  Describing  Funotion  Predictions 

Since  there  are  an  infinite  number  of  nonlinearities  which 
could  be  introduced  into  an  infinite  number  of  feedback  systems, 
the  44  varieties  analyzed  in  this  paper  certainly  do  not  constitute 
a  conclusive  proof  of  the  corrected-conventional  describing  func¬ 
tion's  superior  accuraoy.  But  on  the  basis  of  this  Btudy,  it  cer¬ 
tainly  seems  permissible  to  postulate  that  the  corrected-conventional 
describing  function  can  usually  give  the  control  engineer  a  more 
quantitative  answer  to  his  stability  problem  than  he  vas  able  to 
previously  obtain  from  the  conventional  describing  function,  except 
perhaps  in  a  very  complex  system  vhich  provides  almost  perfect  low 
pass  filtering.  The  reason  for  this  can  be  gleaned  to  some  extent 
from  the  prediction  tables  in  Chapter  IV.  A  close  analysis  of  the 
prediction  tables  in  Chapter  IV  shows  that  the  amplitude  accuracy 
of  the  conventional  describing  function  decreases  as  the  limit- 
cycle  amplitude,  X,  increases.  This  is  so  because,  as  Thaler  and 
Pastel  point  out  (Ref  12:l'f2),  the  conventional  describing  func- 
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tion  'becomeE  less  dependable  as  the  signal  level  increases.  In 
addition,  the  amplitude  accuracy  of  the  conventional  describing 
function  decreases  as  the  order  of  the  system  decreases,  since  less 
filtering  of  the  higher  harmonics  occurs.  And,  as  the  tables  show, 
the  amplitude  accuracy  of  the  conventional  describing  function  also 
seems  to  depend  on  system  type,  decreasing  as  the  system  type  is 
increased.  The  corrected-conventional  describing  functions,  on  the 
other  hand,  seem  to  be  relatively  unaffected  by  changes  in  limit- 
cycle  amplitxide,  system  type,  or  system  order  for  the  intermediate 
order  systems  of  the  study.  This  seems  to  come  about  from  the 
inherent  correction  in  the  variation  afforded  by  the  third  harmonic. 
Thus,  by  using  the  corrected-conventional  describing  function,  the 
engineer  can  remove  the  small-signal  limitation  and  relax  somewhat 
the  sine-wave  input  assumption,  thus  providing  a  more  quantitative 
analysis  tool  lor  his  real-life  problems. 

The  amplitude  accuracy  of  tbe  new  rms  describing  function  is 
also  relatively  insensitive  to  limit-cycle  amplitude  and  system 
order  changes  in  the  intermediate-order  system  range.  However,  it 
is  not  as  accurate  as  the  corrected-conventional  describing  func¬ 
tion,  and  it  seems  to  be  quite  sensitive  to  changes  in  system  type. 

For  the  purpose  of  attaining  accurate  amplitude  predictions 
with  the  describing  function  stability  analysis  technique,  it 
appears  that  the  minimum  average  error  and  equivalent  gain  schemes 
are  completely  out  ol  the  question.  This  is  to  be  expected,  how¬ 
ever,  since  Gibson  introduced  the  minimum  average  error  technique 
only  as  an  example,  and  Prince’s  equivalent  gain  technique  is 
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intended  primarily  for  analytical-graphical,  closed-loop  frequency 
response  analysis.  Graphical  methods  for  producing  the  closed- 
loop  frequency  response  from  any  describing  function  scheme  do 
exist,  however,  (fief  12:195-210)  and  the  corrected-conventional 
describing  function  applied  in  this  manner  might  produce  more  accu¬ 
rate  results  than  Prince's  method,  at  least  for  nonlinearities 
which  do  not  produce  a  phase  shift. 

All  five  types  of  describing  functions  presented  in  this  paper 
predict  the  same  frequency  of  oscillation,  except  when  the  nonlin¬ 
earity  produces  a  phase  shift,  as  the  relay  +  hysteresis  and  relay 
+  dead  zone  +  hysteresis  nonlinearities  in  this  study  do.  In  those 
cases  it  seems  that  the  equivalent  gain  describing  function  is  more 
accurate  for  frequency  prediction,  with  the  accuracy  increasing  as 
the  limit-cycle  amplitude  is  increased.  This  increase  in  accuracy 
with  amplitude  may  "be  due  in  part  to  the  same  sort  of  phenomenon 
explained  by  Johnson  for  the  conventional  describing  function  (Ref 
6:175);  that  is,  as  the  limit-cycle  amplitude  is  increased,  the  inter¬ 
section  of  the  negative  inverse  describing  function  locus  and  the 
direct  polar  plot  locus  become  more  orthogonal  (see  Figures  12  and 
16,  pages  94  and  69),  and  it  is  the  nearly  tangential  intersections 
which  give  misleading  results. 

If  the  corrected-conventional  describing  function  is  to  be 
acceptable  to  the  control  system  engineer,  it  must  be  shown  why,  in 
fact,  this  new  technique  should  produce  more  accurate  amplitude  pre¬ 
dictions  than  the  conventional  describing  function,  or  any  of  the 
others. 
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The  Case  i'or  the  Corrected-Conventional  Describing  Function. 

The  only  two  methods  of  generating  describing  functions  which 
have  met  with  much  success  prior  to  this  study,  insofar  as  predic¬ 
ting  accurate  limit-cycle  amplitudes  is  concerned,  are  the  conven¬ 
tional  and  new  rms  techniques  (Ref  3:1321).  In  order  for  the  con¬ 
ventional  technique  to  be  completely  accurate,  it  is  necessary  for 
the  linear  portion  of  the  loop  to  act  as  a  perfect  low-pass  filter. 
However,  for  intermediate  order  systems,  the  linear  portion  of  the 
loop  is  far  from  being  a  perlect  low-pass  filter.  Consequently, 
some  higher  harmonics,  although  attenuated,  are  led  back  to  the 
input  of  the  nonlinearity  along  with  the  fundamental.  For  inter¬ 
mediate  order  systems,  the  conventional  describing  function  will, 
thereiore,  underestimate  the  amplitude,  X,  of  the  signal  returning 
to  the  nonlinearity. 

The  new  rms  describing  function,  on  the  other  hand,  assumes 
that  their  is  no  filtering  of  the  higher  harmonics  by  the  linear 
portion  of  the  loop,  but  rather,  that  the  output  sinusoid  is  just 
reshaped  into  a  fundamental  frequency  sine  wave  of  equal  energy 
(Ref  2:3cl).  Appendix  F  shows  that  the  new  rms  describing  function, 
in  effect,  takes  as  the  amplitude  of  the  equivalent  output  sine  wave 
from  the  nonlinearity,  the  square  root  of  the  sum  of  the  squares  of 
all  the  Fourier  coefficients  in  the  Fourier  expansion  of  the  actual 
output  of  the  nonlinearity.  Obviously,  some  liltering  of  harmonics 
higher  than  the  fundamental  will  result  in  an  intermediate  order 
system,  so  the  new  rms  describing  : unction  will,  therefore,  over- 
cstiaate  the  amplitude,  X,  ol  the  signal  returning  to  the  non- 


85 


QGc/EE/64-16 


linearity. 

What  actually  happens  to  the  higher  harmonics  lies  between  what 
the  conventional  and  new  rms  methods  assume.  The  real  case  for  the 
corrected-conventional  describing  funotion  is,  therefore,  that  it 
considers  only  the  most  prominent  higher  harmonic,  in  addition  to 
the  fundamental,  by  taking  as  the  amplitude  of  the  equivalent  output 
sine  wave,  the  square  root  of  the  sum  of  the  squares  of  the  first 
and  third  Fourier  coefficients  in  the  Fourier  Series  expansion  of  the 
actual  output  of  the  nonlinearity.  The  corrected-conventional  method 
therefore  assumes  that  the  joint  effect  of  all  of  the  higher  harmon¬ 
ics,  each  attenuated  by  a  different  amount  through  the  linear  portion 
of  the  loop,  can  be  adequately  approximated  by  including  the  unatten¬ 
uated  effect  of  the  third  harmonic  only.  Thus,  the  prediction 
afforded  by  the  corrected-conventional  describing  function  will  lie 
between  that  of  the  conventional  and  new  rms  techniques  and  will 
therefore  be  more  realistic  than  either.  'Die  remits  of  the  exper¬ 
imental  portion  of  this  study  support  this  cc  jion  quite  well. 

Recommendations 

Use  of  the  third  harmonic  for  amplitude  prediction  correction 
of  the  conventional  describing  function  by  the  technique  of  the 
corrected-conventional  describing  function  defined  on  page  26  has  met 
with  much  success,  as  attested  by  this  paper.  Perhaps  by  applying 
some  sort  of  correction  to  the  angle,  as  veil  as  to  the  magnitude  of 
the  conventional  describing  function,  the  loss  in  frequency  prediction 
accuracy  for  nonlinearities  which  produce  phase  shifts  could  be 
reclaimed,  or  even  improved.  That  is,  a  suitable  describing  func- 
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tion  for  future  study  might  "be  defined  by  the  relation 


where  A^,  A^,  B^,  and  are  defined  as  on  page  26.  This  defini¬ 
tion  would  give  the  same  describing  function  as  the  corrected- 
oonventional  definition  when  the  nonlinearity  output  is  naturally 
odd  periodio,  hut  would  modify  the  angle  of  the  corrected-conven- 
tional  describing  function  when  the  output  of  the  nonlinearity  is 
not  naturally  odd  periodic.  Although  the  resulting  describing 
functions  might  appear  very  complex,  they  could  be  tabulated  easily 
with  the  aid  of  a  digital  computer  and  plotted.  Once  an  accurate 
plot  is  obtained,  it  makes  little  difference  how  complex  the  rela¬ 
tionship  was  from  which  it  came. 

The  success  of  the  correoted-conventional  describing  function 
in  improving  the  prediction  of  the  amplitude  of  self-sustained 
oscillations  in  a  nonlinear  feedback  system  perhaps  warrants  further 
exploration.  The  corrected-conventional  describing  functions  of 
nonlinearities  not  included  in  this  paper  should  be  derived,  cata¬ 
logued,  and  tested  for  accuracy.  Such  additional  nonlinearities 
might  include  negative  defficiency  (Ref  12:150),  backlash,  nonlin¬ 
earities  described  by  algebraic  equations  such  as  y  «  x-|xj  or 
y  -  x^,  and  nonlinearities  described  by  differential  equations. 

Finally,  it  might  be  interesting  to  see  which  type  of  describ¬ 
ing  function  can  most  accurately  predict  the  closed  loop  frequency 
response  and  jump  resonance. 
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All  of  the  questions  posed  "by  the  preceeding  recommendations 
must  be  answered  before  the  question  "Which  describing  function  is 
best?"  can  be  fully  answered.  However,  based  on  the  results  of 
this  paper,  it  seems  safe  to  recommend  the  following  procedure  for 
treating  the  stability  analysis  of  intermediate-order  systems: 

(1)  For  dead  zone,  dead  zone  with  saturation,  ideal  relay,  and 
relay  with  dead  zone  nonlinearities,  use  the  corrected-conventional 
describing  function  for  the  best  results  in  both  amplitude  and 
frequency  prediction  of  limit  cyoles. 

(2)  For  relay  with  hysteresis  and  relay  with  both  hysteresis 
and  dead  zone,  use  the  oorreoted-conventional  describing  function 
for  the  amplitude  prediction  and  the  equivalent  gain  describing 
function  for  the  frequency  prediction. 

(3)  For  saturation  nonlinearities,  use  the  conventional  de¬ 
scribing  function  for  small-signal  analysis  and  the  corrected- 
conventional  describing  function  for  large-signal  analysis;  although 
not  muoh  will  be  lost  by  using  the  corrected-conventional  describing 
function  exclusively. 
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Appendix  A 

Derivations  of  the  Correot6d-Conventlonal 
Describing  Functions  Used  in  this  Study 


The  definition  of  the  correoted-conventional  describing  function 
given  in  Eq  (16)  on  page  26  is  quite  general,  but  the  definition  can 
often  be  simplified  by  taking:  advantage  of  certain  symmetries  in  the 
output  waveform.  For  odd-periodic,  odd-harmonic  outputs  from  the 
nonlinearity,  Eq  (16)  becomes  simply  (nondimensionalised) 


(T2+  T2)^ 
V  V  1  3  ; 

K  “  KX 


(18) 


o  p  JL 

where  «  Fundamental  amplitude  «  ^  )7  «  Bx,  since  -  0* 

2 

-  3rd.  harmonic  amplitude  «  (A£  •+  B 
X  «  Amplitude  of  the  sinusoidal  input  tc  the  nonlinearity; 

K  ■  Gain  of  the  nonlinearity. 

Obviously,  the  derivation  of  the  corrected-conventional  describing 
function  is  similar  to  that  of  the  conventional,  exoept  that  the 
amplitude  of  the  third  Fourier  harr.onic  must  be  determined  as  well 
as  the  fundamental.  As  a  first  example,  oonsider  the  case  of 
saturation  combined  with  dead  zone. 


B 


3’ 


since 


Saturation  Combined  with  Dead  Zone 

As  Figure  19  on  the  following  page  shows,  the  output  of  the  non¬ 
linearity  is  both  odd  periodic  and  odd  harmonic;  therefore,  the  coef¬ 
ficients  of  the  cosine  terms  in  the  Fourier  expansion  of  the  output 
waveform  are  all  zero,  as  are  the  coefficients  of  all  the  even  sine 
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zone  is  therefore 


F 

K 


1 

TT 


(2ch-b  -  2*^  +  sin  2^  -  sin  2oc^)2  +  (-|  sin  2 ctg  - 
sin  40^  +  sin  4«d  -  y  sin^.  cos  ocg  + 


4  sin  2a, 
3  d 


-l  £  •J  p 

y  sin  <xd  cos  a.^) 


£ 


(23) 


The  (nondimensional)  corrected-conventional  describing  function 
for  saturation  can  he  found  by  letting  oC^  ■  0  in  Eq  (23);  that  is, 


F 

K 


+  sin  2c*b) 


2 


+ 


(*|  sin  2^  -  sin 

0  2 
-j  sin  2oc-b  sin  ocb) 


(24) 


Similarly,  by  letting  oia  «TT/2  in  Eq.  (23)  ,  one  can  find  the  (nondi- 
mensional)  corrected-conventional  describing  function  for  dead  zone 
to  be 


F 

K 


1 

ir 


pc  *1 

(If-  2oCd  -  sin  2«d)  +  (“  3  sin  2°^  +  -p  Bin 


8  • 

-r  s xu  2a, 
3  d 


.  2^  x2 

sxn  0^) 


4cSi  + 

* 


(25) 


Eqe  (24)  and  (25)  can  also  be  obtained  by  a  Fourier  analysis  of  the 
output  wave  of  the  nonlinearity  as  was  done  for  saturation  combined 
with  dead  zone. 


Relay  with  Dead  Zone  and  Hysteresis 

Using  the  same  technique  as  in  the  preceding  derivation,  one 
can  derive  the  (nondimensional)  corrected-conventional  describing 
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function  for  a  relay  with  dead  zone  and  hysteresis,  and  then  simplify 
the  result  for  the  specific  cases  of  a  relay  with  dead  zone  only,  a 
relay  with  hysteresis  only,  and  an  ideal  relay. 

Figure  20  on  the  following  page  shows  that  the  output  of  the 
relay-dead  zone-hysteresis  nonlinearity  is  not  symmetrical  about  the 
origin.  Eowever,  odd-periodic  and  odd-harmonic  symmetry  can  he 

recovered  by  a  phase  shift  of  the  axis  by  the  amount  (7T/2 - — - ). 

This  angle  is  therefore  the  angle  of  any  describing  function  which 
is  generated  by  a  scheme  which  does  not  introduce  a  phase  shift  when 
applied  to  an  odd-periodic,  odd-harmonic  function.  Since  all  of  the 
describing-function  generating  schemes  studied  in  this  paper  meet 
this  requirement,  the  expression  for  the  angle  of  the  corrected- 
conventional  describing  function  for  a  relay  with  dead  zone  and 
hysteresis  is 


K 


TT  Ct+/g 
2  2 


(26) 


After  the  axis  is  shifted,  the  output  waveform  is  odd-periodic 
and  odd-harmonic.  Therefore,  all  of  the  cosine  coefficients  and  the 
even  sine  coefficients  are  zero,  and  the  amplitude  of  the  funda¬ 
mental  (integrating  with  the  axis  shifted  to  -  /^S)  is 


f? 

Yi  “  4  JQ  3^*) 


T r 

K  sin G  d  9 


(27 


or 
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From  Eqa  (18),  (2'9),  and  (31)  the  magnitude  of  the  (nondimen- 
sional)  corrected-conventional  describing  function  for  a  relay  vith 
dead  zone  and  hysteresis  is 


or,  substituting  X 


d+h 

2  sin*  ’ 
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The  (nondimensional)  corrected-oonventional  describing  function 


for  a  relay  with  dead  zone  only  can  be  determined  by  letting  h  -  0 


in  the  foregoing  analysis}  then, 


ot  -  sin”1 

cJl.  CL 


/5  -  tt-  Sin”1  ^  .  7T-  *d 


So  Eq  (26)  becomes 


rr  <*d 


zl-f- 


+  7 7  -  a.. 


Therefore  the  describing  function  for  a  relay  with  dead  zone  is 


determined  from  Eq  (33)  to  be 


Hd  4  . 

-rr-  -  -?  83J 


K  7 T 


i  *  -5(1  - 


4  sin  «c.  )  ‘ 

d 


Similarly,  by  letting  d  -  0 


a.  six,-1  ^  =  «h 


/3  -  7T-  sin”1  ~  -  7T-f  a:h 


And,  by  substituting  these  values  of  et  and  /S  into  Eqs  (26)  and  (32), 


one  can  find  the  (nondimensional)  corrected-conventional  describing 


function  for  a  relay  with  hysteresis  to  be 


d-Zi-A 


k|_  A.  (i°)4 
K  I  to  v  9  ' 


■ 
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Furthermore,  for  an  ideal  relay  OL^  -  0,  and  therefore 

!-*<¥>*  <«> 

Eqs  (37) »  (40) ,  (4l)»  and  (42)  can  also  be  found  by  a  Fourier 
analysis  of  the  output  vavef orm  of  the  respective  nonlinearities  as 
•was  done  for  saturation  combined  with  dead  zone  and  the  relay  with 
dead  zone  and  hysteresis. 

Plots  of  all  of  the  describing  functions  derived  in  this  appen¬ 
dix  can  be  found  in  Appendix  D. 
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Appendix  B 


Derivations  of  the  Hew  HMS 


Functions  which  Introduce 


As  pointed  out  on  pages  6  and  8  of  the  introduction,  one  way 


of  extending  the  new  rms  describing  function  to  polar-symmetrical 


nonlinearities  which  are  not  Bingle-valued  is  to  let  the  angle  of 


the  describing  function  be  the  angle  of  phase  shift  necessary  to 


make  the  output  of  the  nonlinearity  odd  periodic.  This,  incidently, 


is  the  same  angle  that  is  obtained  for  the  conventional  describing 


function  through  a  Fourier  analysis.  Again,  as  with  the  corrected- 


conventional  describing  function  derivations  in  Appendix  A,  one  can 


take  advantage  of  the  quarter-wave  symmetry  which  occurs  after  the 


shift  in  axes  to  modify  Eq  (15)  on  page  25  for  simpler  integration. 


For  example,  after  shifting  the  axes  of  the  output  of  the  relay  + 


dead  zone  +  hysteresis  nonlinearity  shown  in  Figure  20  on  page  96 


to  the  output  becomes  odd  periodic,  and  the  angle  of  the  new 


rms  describing  function  is,  therefore, 


Then,  the  magnitude  of  the  new  rms  describing  function  is 


\  X2  sin2©  d© 

Jo 


K2  d © 


I 
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because  of  the  quarter  wave  symmetry. 


x 

K2® 

2 

n  +°L=£ 

2  2 

X2  (f  -  iisJfi.) 

IT 

2 

0_ 

Integrating,  one  will  obtain 


(45) 


So  the  (nondimensional)  new  rms  describing  function,  with  angle 
expressed  in  exponential  form,  is 


N 

K 


C< 

2  ' 


(46) 


or,  substituting  X 


dth 
2  sinct* 


N(d+h) 

K 


2 

TT 


{$-0 0  sin  * 


*  dO 


T_«±A 


(47) 


Ey  letting  d  »  0  in  Figure  20,  one  can  find  the  (nondimensional) 
new  rms  describing  function  for  a  relay  with  hysteresis.  Then 


.  -1  h  _  „ 

«  *  61n  2X  *  *h 

/S  =  7T  -  sin-1  ~  =  TT  +  if 


N  V2-  "A 
K  =  T  6 


(33) 

(39) 


(46) 
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Appendix  C 

Experimental  Procedure  and  Analog  Computer  CircuitB 


In  order  to  add  credibility  to  the  experimental  results,  as 
well  as  to  allow  the  reader  to  repeat  the  experiment,  the  experi¬ 
mental  procedure  used  by  the  author  is  explained  in  the  following 
discussion. 


Procedure 

Figure  21  on  the  following  page  shows  the  computer,  plotter, 
and  test  equipment  arrangement  used  in  the  experimental  determina¬ 
tion  of  the  amplitude  and  frequency  of  the  self-sustained  oscilla¬ 
tions  in  the  nonlinear  feedback  systems  of  the  study.  Figures  23 
through  29  on  pages  107  through  111  show  general  analog  computer 
circuits  equivalent  to  the  block  diagram  arrangement  in  Figure  22 
for  each  of  the  linear  forward  transfer  functions,  G^s)  through 
G^(s).  The  general  circuit  for  G^s)  is  the  sane  as  that  of  G^(s) 
in  Figure  25  >  but  with  the  outer  loop  gain  increased  by  a  factor 
of  3.475*  Similarly,  the  general  circuit  for  G71(s)  is  the  same 
as  that  of  G^(s)  in  Figure  29  >  but  with  the  outer  loop  gain  in¬ 
creased  by  3-333.  Figures  30  through  36  on  pages  111  through  115 
contain  the  inferential  relay  mechanizations  of  the  various  non- 
linearities  studied.  These  inferential  relay  circuits  were  intro¬ 
duce!  m  place  01  the  box  labled  "F'  in  each  of  the  general  cir¬ 
cuits  lor  tr.e  specific  tests  required.  The  procedure  for  measure¬ 
ment  of  the  resulting  oscillations  was  as  follows  (See  Figure  21): 
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Figure  21 

Analog  computer  test  arrangement  for  the  measurement  of  the  ampli¬ 
tude  and  frequency  of  self-sustained  oscillations:  (l)  REAC  Elec¬ 
tronic  Analogue  Computer,  Kodel  C-101,  Reeves  Instrument  Corpora¬ 
tion;  (2)  Differential  Relay  Unit;  (3)  Yariplotter  Plotting  Board, 
Kodel  1100— E,  Electronic  Associates,  Inc.;  (4)  Dynagraph  Recorder, 
Type  KC,  Offner  Electronics,  Inc.;  (5)  Digital  Voltmeter,  Kodel 
4011-4,  Beckman  Industries,  Inc. 


G0C/W64-16 


(1)  The  diiferential  relay  circuit  and  general  circuit  lor  a 
particular  case  to  be  tested  were  wired  together  on  the  READ  patch- 
panel,  alter  proper  amplitude  scaling  based  on  the  predicted  ampli¬ 
tude.  (Note  that  the  figures  in  this  appendix  have  been  labled  in 
terms  of  k,  the  amplitude  scale  2 actor,  in  order  to  facilitate 
scaling  each  run.)  The  potentiometers  on  the  RLAC  were  then  set  to 

f 

four-place  accuracy  with  the  digital  voltmeter. 

(2)  The  R&AC  was  switched  to  "Operate"  and  the  ensuing  oscil¬ 
lations  were  recorded  on  both  the  Variplotter  and  the  Dynagraph. 

An  initial  value  of  c  which  was  greater  than  the  dead  zone  was 
required  to  initiate  oscillations  in  the  systems  that  had  nonlin¬ 
earities  with  dead  zone. 

(3)  After  the  oscillation  had  stabilized  for  a  few  cycles,  the 
RLAC  was  switched  back  to  "Reset." 

(4)  The  amplitude  of  the  oscillation  was  then  accurately  deter¬ 
mined  by  measuring  the  positive  and  negative  initial  condition  vol¬ 
tage  (on  the  integrator  from  which  X  was  being  recorded)  necessary 
to  move  the  pen  oi  the  Variplotter  precisely  to  the  level  0.  the 
positive  and  negative  peaks  of  the  previously  recorded  seli-sus- 
tained  oscillation.  This  initial  condition  voltage  was  measured 

to  four-place  accuracy  with  the  digital  voltmeter.  The  amplitude 
of  the  oscillation  was  then  taken  as  the  avenge  of  the  magnitude 
of  the  positive  end  negative  peak  so  measured.  By  using  this 
technique,  one  can  eliminate  any  variation  in  the  linearity  of 
the  Variplotter,  and  very  accurate  results  can  be  obtained. 
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(5)  The  frequency  of  the  oscillation  was  accurately  deter¬ 
mined  by  measuring  the  period  of  a  previously  recorded  self-sus¬ 
tained  oscillation  from  a  Dynagraph  recording.  Since  the  Qyna- 
graph  was  run  at  one  of  its  maximum  speeds  during  the  test  (100 
or  250  mm/sec)  and  the  frequency  of  oscillation  was  quite  low  (0.5 
to  3.3  radians/sec) ,  the  measurement  of  the  period  of  the  oscil¬ 
lation  was  also  quite  precise.  The  frequency  in  radians/sec  is 
then  2TTdivided  by  the  period. 

Verification 

In  order  to  verify  the  accuracy  of  this  experimental  procedure, 
one  of  the  runs,  G^(s)  plus  saturation  and  dead  zone,  was  programmed 
and  run  on  the  digital  computer.  G^(s)  was  chosen  as  the  linear 
function  to  be  checked  because  it  was  the  worst  of  all  of  the  analog 
circuit  designs  due  to  the  high  loop  gains.  The  dead  zone  and 
saturation  nonlinearity  was  chosen  because  it  was  one  of  the  more 
complex  nonlinearities  studied  from  a  standpoint  of  simulation 
difficulty,  and  it  was  felt  that  the  sensitivity  of  the  differen¬ 
tial  relays  might  introduce  a  small  error  in  the  simulations,  and 
this  error  would  be  most  prominent  in  a  nonlinearity  with  several 
relays. 

A  numerical  integration  time  increment  of  0.004  seconds  was 
used  for  the  digital  computer  study.  The  results  indicated  a  limit 
cycle  at  the  output  with  an  amplitude  of  12.632,  and  a  frequency 
of  1.622  rad/sec.  The  analog  computer  run  gave  an  amplitude  of 
12.66  and  a  frequency  of  1.616  rad/sec  for  the  sane  problem. 


106 


CCC/EE/64-16 


Since  the  iiiierence  between  the  digital  and  analog  results 
was  less  than  half  a  percent  (see  Table  XXVIII  below)  in  what  is 
considered  to  be  one  of  the  least  accurate  analog  runs,  it  is  felt 
that  the  validity  of  the  entire  set  of  analog  computer  results  is 
sufficient  for  a  basis  of  comparison  of  the  various  describing 
functions. 


Table  XXV.III 

Comparison  of  Analog  and  Digital  Simulations 


Amplitude  Frequency 


Digital  Computer  Result 

12.832 

1.822 

Analog  Computer  Result 

12.88 

1.816 

‘/e  Difference  in  Results 

0.3741 

0.3293 
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Figure  24 

General  analog  computer  circuit  for  the  determination  of  the 
amplitude  and  frequency  of  sustained  oscillations  in  a  dosed 
loop  vith  linear  system  02(a)  and  a  nonlinearity. 

Transfer  function)  °2(«)  -  .(W)(s+2)  -  ffj} 

y  o 

Mechanized  differential  equation!  Do  -  l8y  -  3D  o  -  2Dc 
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Figure  26 

General  analog  computer  circuit  for  the  determination  of  the 
amplitude  and  frequency  of  sustained  oscillations  in  a  closed 
loop  with  linear  system  Q^(s)  and  a  nonlinearity. 

Transfer  function,  O^a)  . 

Mechanised  differential  equation! 

D50  -  46 50* 3 J  -  15D*c  -  85D3o  -  225D2C  -  274D0  -  120c 
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Figure  27 

General  analog  conputer  circuit  for  the  determination  of  the 
amplitude  nn*  frequency  of  sustained  oscillations  in  a  closed 
loop  with  linear  system  Q,-(s)  and  a  nonlinearity. 


Transfer  function,  G^e)  m 

Xechanized  differential  equation: 

D5c  -  301y  -  1QD4o  -  35D3c  -  50D2c  -  24Dc 


*7  #12 


->k6  — 


'r\ v 

(.6  2  5? — ^20 


Figure  25 

General  analog  computer  circuit  for  the  determination  of  the 
amplitude  and  frequency  of  sustained  oscillations  din  a  closed 
loop  with  linear  system  G^(»)  and  a  nonline  arity. 

Transfer  function:  G.-(b)  -  (a+o.l)  -  uv 

6  ^ 

Kechanized  differential  equation: 

D5c  -  242.97  +  24.29y/D  -  12D4o  -  461^0  -  60D2c  -  25Dc 
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Figure  29 

General  analog  computer  circuit  for  the  determination  of  the 
amplitude  and  frequency  of  sustained  oscillations  in  a  closed 
loop  with  linear  ayaten  G^(a)  and  a  nonlinearity. 

OiQ  £.**£.  /  —  _  \  2  —  /  \ 


Transfer  function,  0,(s)  , 

7  a3(a+10)2 

Mechanized  differential  equation, 

D3o  -  248.626y  ♦  497-252  +  248.626  -  20D2c  -  lOODc 

D2 


y 

kx 

DH3I  A 


i  -100 


rigure  30  j 

i 

) 

Mechanization  of  dead  zone  (d  -  6,  I  -  0.5)  with  differential 
relays:  (a)  transfer  characteristic,  (h)  analog  cocputer  circuit-! 

i 

(FroE  Ref  5  : 192)  ! 
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Figure  31 

Mechanisation  of  saturation  (S  -  3,  K  -  0.5)  with  differential 
relay* x  (a)  transfer  characteristic,  (h)  analog  computer  circuit. 

(Fro*  Ref  5  ,188) 


slope  /I 


figure  32 

Mechanization  of  saturation  with  dead,  zone  (d  -  2,  S  -  4,  K  -  l)  with 
differential  relays:  (a)  transfer  characteristic,  (b)  analog  computer 
circuit. 
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Figure  33 


Mechanisation  of  an  idaal  relay  (t  -  1)  differential  relays: 

(a)  transfer  characteristic,  (h)  analog  oonputor  circuit. 

(Fro*  Eef  5  :190) 


Figure  34 

Mechanization  of  a  relay  with  hysteresis  (h  -  6,  K  -  l)  with 
differential  relays:  (a)  transfer  oh&raoteristic,  (h)  analog 
conputer  circuit.  - - - - - — . — - 
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figure  35 

Mechanisation  of  a  relay  vith  dead  zone  (d  «  3,  K  -  l)  vith 
differential  relays:  (a)  transfer  characteristic,  (b)  analog 
computer  circuit. 
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Appendix  D 

Plots  of  the  Corrected-Conventional  describing  Functions 

Figures  37  through  44  on  the  following  pages  are  plots  of  the 
corrected-conventional  describing  functions  derived  in  Appendix  A 
and  listed  in  Tables  II  through  VIII  in  Chapter  III.  The  plots 
are  normalized  as  were  those  in  Chapter  IV,  but  the  ordinate  and 
the  abscissa  are  the  reciprocal  of  those  in  the  graphs  of  Chapter 
IV.  This  is  done  so  that  the  entire  graph  can  be  plotted,  without 
having  one  scale  or  the  other  become  infinite.  (The  ideal  relay 
plot  is  an  exception,,  but  it  is  merely  a  straight  line  relation¬ 
ship.)  Therefore  each  curve  is  a  complete  presentation  of  the 
corrected-conventional  describing  function  for  that  nonlinearity. 
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Figure  41 


Angle  of  the  corrected-conventional  describing  function  of  a 
relay  vith  hysteresis,  (the  magnitude  variation  oi  the  describ¬ 
ing  function  of  a  relay  with  hysteresis  is  the  same  as  for  an 
ideal  relay,  and  therefore  the  curve  on  the  preceding  page  can 
he  used)  where 
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Figure  43 

A  family  of  curves  for  the  magnitude  of  the  non-dimenaionalized 
corrected-conventional  describing  function  of  a  relay  vith  dead 
zone  and  hysteresis,  where 
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Appendix  E 

Other  Describing  Function  Derivations 

Although  each  of  the  35  describing  functions  listed  in  Tables 
II  through  VIII  in  Chapter  III  vas  derived  by  the  author,  only  the 
derivations  of  the  important  now  describing  functions  are  included 
in  this  paper  (Appendices  A  and  B) .  The  following  tables  are  in¬ 
cluded,  therefore,  for  the  scrutinizing  reader  who  desires  to  check 
the  derivations  of  the  other  describing  functions. 


Table  XXIX 

References  where  Derivations  of  the  Conventional 
Describing  Functions  Used  in  this  Study  Kay  Be  Pound. 


Nonlinearity 

Reference 

Dead  Zone 

1:434 

Saturation 

1:435 

Deal  Zone  +  Saturation 

1:436 

Ideal  Relay 

9:457 

Relay  +  HysteresiE 

11:740 

Relay  +  Dead  Zone 

9:458 

Relay  +  Dead  Zone  +  hysteresis 

9:459 
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Table  XXX 

References  vhere  Some  of  the  Unconventional  Describing 


Functions 

Used  in  this  Study  Fay 

Be  round. 

Describing  Function 

Hbnlinearity 

Reference 

Mew  RMS 

Saturation 

3: 1321 

Ideal  Relay 

2:332 

Relay  +  Dead  Zone 

2:382 

Minimum  Average  Error 

Ideal  Relay 

2:381 
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Appendix  F 

Analysis  of  the  New  HKb  Describing  Function 


In  Chapter  III  the  new  nns  describing  function  is  defined  as 


U5) 


where  X  -  the  amplitude  of  the  sinusoidal  input  to  the  nonlinearity 
and  y(t)  -  the  actual  output  from  the  nonlinearity.  By  comparison 
with  Eq  (l)  in  Chapter  I,  Eq  (15)  indicates  that  the  amplitude  of  the 
equivalent  output  sine  wave  from  the  nonlinearity  has  been  chosen  as 


(49) 


That  is,  the  equivalent  output  sine  vave  has  the  same  rms  value  as 
the  actual  output  of  the  nonlinearity.  If  y(t)  is  an  odd-harmonic, 
odd-periodic  function,  as  is  often  the  case,  and  is  expanded  in  a 
Fourier  Series,  then  Eq  (49)  becomes 

~1  f2Tr  2  1  ^ 

I  -  —  V  (y]L  sin (jiX  +  sin  3wt  +  sin  5*ot  +  ....)  d^ 

1  f2^  ?  2 

-  —  \  (Yn  sin  cot  +  2Y..Y  '  sinuJt  sin  3^t  +  2Y.YC  oinajt  sin  5 «4t 
"n  J0  1  id  1  5 

,2  .2. 


+  ...  +  Y^  sin  i'-'t  +  2Y^Y^.  sin  3e2t  sin  cjcJt  +. 

£ 


+  Y^2  cin2^*  *■  ••••)  d<*Jtj 


'2tt 


but  \  sin  owt  sin  ncot  durt  «  0 

Jo 


if  B  / 


n 


(50) 

(51) 
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(FT  2 

and.  \  sin  nwt  dcot  *•  tt 

Jo 


Therefore,  |t|  .  [|  (T**  ♦  ^  + 


(52) 


(53) 


Stating  the  above  analysis  in  vords,  the  nev  rms  describing  function 
in  effect,  takes  as  the  amplitude  of  the  equivalent  output  sine  vave 
from  the  nonlinearity,  the  square  root  of  the  sum  of  the  squares  of 
an  the  Fourier  coefficients  in  the  Fourier  Series  expansion  of  the 
actual  output  of  the  nonlinearity. 


